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ABSTRACT 



Using the largest high-resolution spectroscopic sample to date of young, very 



, low mass stars and brown dwarfs, we investigate disk accretion in objects rang- 

r ^ \ ing from just above the hydrogen-burning limit all the way to nearly planetary 

masses. Our 82 targets span spectral types from M5 to M9.5, or masses from 



0.15 M down to about 15 Jupiters. They are confirmed members of the p Ophi- 
uchus, Taurus, Chamaeleon I, IC 348, R Coronae Australis, Upper Scorpius and 
TW Hydrae star-forming regions and young clusters, with ages from <1 to ~10 
Myr. The sample contains 41 brown dwarfs (spectral types >M6.5). We have 
previously presented high-resolution optical spectra for roughly half the sample; 
i/-) | the rest are new. This is a close to complete survey of all confirmed brown dwarfs 

known so far in the regions examined, except in p Oph and IC 348 (where we are 
limited by a combination of extinction and distance). We find that: (1) classical 
T Tauri-like disk-accretion persists in the sub-stellar domain down to nearly the 



O 

deuterium-burning limit; (2) while an Ha 10% width > 200 km s 1 is our prime 
accretion diagnostic (following our previous work), permitted emission lines of 
Ca II , O I and He I are also good accretion indicators, just as in CTTs (we cau- 
tion against a blind use of Ha width alone, since inclination and rotation effects 



on the line are especially important at the low accretion rates in these objects); 
(3) the Ca II 8662A line flux is an excellent quantitative measure of the accretion 
rate in very low- mass stars and brown dwarfs (as in higher- mass CTTs), corre- 
lating remarkably well with the M obtained from veiling and Ha-modeling; (4) 
the accretion rate diminishes rapidly with mass - our measurements support pre- 
vious suggestions that M oc M. 2 (albeit with considerable scatter), and extend 
this correlation to the entire range of sub-stellar masses; (5) the fraction of very 
low-mass stellar and brown dwarf accretors decreases substantially with age, as 
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in higher-mass stars; (6) at any given age, the fraction of very low- mass stellar 
and substellar accretors is comparable to the accretor fraction in higher-mass 
stars; and (7) a number of our sources with infrared excesses arising from dusty 
disks do not evince measurable accretion signatures, with the incidence of such a 
mismatch increasing with age: this implies that disks in the low mass regime can 
persist beyond the main accretion phase, and parallels the transition from the 
classical to post-T Tauri stage in more massive stars. These strong similarities at 
young ages, between higher-mass stars and low-mass bodies close to and below 
the hydrogen-burning limit, are consistent with a common formation mechanism 
in the two mass regimes. 

Subject headings: stars: low-mass, brown dwarfs - stars: pre-main sequence - 
stars: formation - circumstellar matter - planetary systems - techniques: spec- 
troscopic 



1. Introduction 

There has been much recent interest and activity in investigating the origin and early 
evolution of sub-stellar bodies. Some theorists, most recently Padoan & Nordlund (2004), 
have suggested that objects all the way from solar-mass stars to brown dwarfs form the same 
way, via 'turbulent fragmentation'. Others have proposed that brown dwarfs are in fact 
"stellar embryos", ejected from newborn multiple systems through dynamical gravitational 
interactions before accreting sufficient mass to become full-fledged stars (Reipurth & Clarke 
2001; Bate et al. 2003). A detailed comparison between the properties of brown dwarfs and 
stars, in their infancy, can help distinguish between these scenarios. 

Recent observations make a compelling case that objects near and below the sub-stellar 
boundary undergo a T Tauri phase, remarkably similar to their higher mass counterparts. 
For example, excess emission in the near- and mid-infrared, indicative of dusty disks, appears 
to be common in young brown dwarfs (e.g., Muench et al. 2001; Natta et al. 2002; Jayaward- 
hana et al. 2003). So do spectroscopic signatures of disk-accretion, such as broad, asym- 
metric Ha lines (Jayawardhana, Mohanty & Basri 2002, 2003 - hereafter JMB02, JMB03; 
Muzerolle et al. 2003 - hereafter MHCBH03; Natta et al. 2004 - hereafter N04). Irregular 
high-amplitude photometric variations, seen in some brown dwarfs, may also be related to 
accretion (Scholz & Eisloffel 2004). What's more, a number of objects around the hydrogen 
burning limit exhibit forbidden line emission, usually associated with jets and winds in T 
Tauri and Class I sources (e.g., Fernandez & Comeron 2001; Barrado y Navascues, Mohanty 
& Jayawardhana 2004; MHCBH03; N04; Barrado y Navascues & Jayawardhana 2004). 
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In this paper, we greatly extend the investigation of accretion in the very low mass 
stellar and sub-stellar regimes, with the largest sample of high-resolution optical spectra for 
such objects to date. Our sample includes more than 80 targets, extending from just above 
the hydrogen burning limit all the way down to masses approaching the planetary domain. 
Thus, we are able to explore accretion in very low mass stars and brown dwarfs over an order 
of magnitude in mass and several orders in derived accretion rates, look for correlations with 
parameters such as age and disk-related infrared excess, and compare these properties to 
those of higher mass T Tauri stars in order to probe the formation mechanism. 

Before proceeding further, it is important to emphasize two points that bear on the 
question of brown dwarf origins. The first concerns the inherent uncertainty in certifying sub- 
stellar status in young objects. In the field, the presence of Lithium is a litmus test for brown 
dwarfs, with objects evincing Li I in their spectra bound to be less massive than ~0.06 M Q 
, and thus indubitably sub-stellar. The test fails at ages of a few Myrs, however, since even 
stars have not yet depleted their Lithium content over such short timescales. In the absence 
of such a direct test, sub-stellarity in young objects is inferred via techniques that are highly 
dependent on theoretical models. At present, the latter models are quite poorly anchored 
in the young very low-mass stellar and substellar regimes, given the paucity of empirically 
determined parameters (temperature, mass, radius) for such objects. Consequently, the 
masses determined from these models are subject to significant uncertainties; combined 
with reasonable observational errors, this can lead to mass offsets of order 50%. Under the 
circumstances, at early ages, no current technique can distinguish between brown dwarfs 
and very low-mass stars that lie close the stellar/sub-stellar boundary: a young "star" with 
inferred mass 0.1 M may well be a 0.07 M brown dwarf instead, and vice versa. For 
ease of nomenclature, one may certainly assume a working definition of sub-stellarity: e.g., 
within the context of the widely-used Lyon evolutionary models, a spectral type of ~M6 is 
usually adopted as the dividing line between stars and brown dwarfs at a few Myrs (here 
we conservatively adopt M6.5; see further below). It must be kept in mind, though, that 
any such boundary is only indicative; the true division may lie, for instance, anywhere from 
M5.5-6.5. Debates over whether an individual young object, at say M6, is "really" a star or 
a brown dwarf are currently impossible to settle, and thus meaningless. 

Given this difficulty in distinguishing between the stellar and sub-stellar domains at 
early ages, how does one choose a sample to explore the formation mechanism of brown 
dwarfs? The first solution is to include a large number of relatively late-type sources that 
are likely to be sub-stellar in spite of the above uncertainties: e.g., a few Myr-old M8 object, 
with inferred mass ~0.03 M Q , is quite certain to be a brown dwarf even if its true mass is a 
factor of 2 higher. Our present sample satisfies this criterion. The second, more important 
answer is that a distinction between very low mass stars and brown dwarfs is not truly 
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relevant to the formation question. The 'ejection' and 'turbulent fragmentation' hypotheses 
seek to explain how very low mass objects can form out of molecular clouds, when the average 
thermal Jeans mass in the latter is usually of order a solar-mass. This is a problem, however, 
that is not restricted solely to the formation of brown dwarfs, but extends to very low mass 
stars as well. If we do not know the formation mechanism of a 0.08 M brown dwarf, we are 
also in the dark about the origins of a 0.1 M star. For example, in the simulations by Bate 
et al. (2003), brown dwarfs and very low mass stars are produced through ejection. In this 
context, the statement often made that brown dwarfs form (or do not form, depending on 
one's preferred hypothesis) like stars draws a spurious distinction at the hydrogen-burning 
limit: no formation scenario put forward so far depends physically on whether an object will 
eventually fuse hydrogen or not. The real physical distinction lies between higher-mass solar- 
type stars, whose formation appears well addressed by the thermal Jeans mass formulation, 
and bodies - both stars and brown dwarfs - with masses that are lower by an order of 
magnitude (or more). It is the formation of all the latter objects that is the fundamental 
question. To probe this issue, it is acceptable to lump together young objects that are both 
close to and below the stellar /sub-stellar boundary, without having to precisely determine 
sub-stellar status; the crucial comparison is between this group and significantly higher-mass 
T Tauri stars, in order to explore whether the same formation process is at play in both mass 
regimes. We carry out such comparisons in the present work. 

In view of the above discussion, we adopt the following nomenclature in this paper. 
"Very low mass stars" (VLMS) are defined as objects with spectral types in the range M5 
- <M6.5. At the 1-10 Myr ages we are concerned with here, this formally corresponds to 
masses, from Lyon theoretical tracks, of 0.15 > .M* > 0.075 M Q (Baraffe et al. 1998). "Brown 
dwarfs" (BDs) are defined as young sources with spectral types >M6.5, corresponding to 
theoretical masses < 0.075 M Q . "Low-mass objects" will refer to the combined VLMS and 
BD sample, i.e., all young sources that are M5 or later (< 0.15 M ). Finally, the latter 
sample will often be compared to "higher-mass stars" , by which we mean young stars with 
spectral types ~ K0 - M4, corresponding to theoretical masses 2 > M* > 0.25 M at a few 
Myrs. The low and high mass ranges we choose allow clear comparisons between stars close 
to a solar mass, and objects smaller than a solar mass by an order of magnitude or more. 



2. Sample and Observations 

Our 82 targets are located in the p Oph, Taurus, IC348, Cha I, R Coronae Australis (R 
CrA), Upper Sco and TW Hydrae star- forming regions and young clusters. These regions 
span ages from <1 Myr to ~10 Myr; the objects range in mass from about 0.15 M down to 
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~15 Mj (age and mass determinations discussed in §3). The sources are all bona-fide cluster 
members: they are selected from low-resolution spectroscopic studies which have already 
established their membership in the respective regions (Luhman & Rieke 1999; Briceno et 
al. 2002; Luhman et al. 2003a, b; Comeron, Neuhauser & Kaas 2000; Fernandez & Comeron 
2001; Barrado y Navascues, Mohanty & Jayawardhana 2004; Gizis 2002; Mohanty et al. 
2004a; Martin, Delfosse & Guieu 2004). We further support youth and membership through 
the detection of lithium (Li I A6708; except in the coolest /faintest objects, which have very 
low S/N in the blue), radial velocities, and the presence of narrow atomic alkali lines (Na I , 
K I ) indicative of low gravity. We have previously presented high-resolution optical spectra 
for 33 of these targets - 4 in p Oph, 1 in R CrA, 3 in the TW Hydrae Association (TW HyA), 
4 in Taurus, 10 in IC348 and 11 in Upper Sco (JMB02; JMB03; Mohanty, Jayawardhana & 
Barrado y Navascues 2003; Barrado y Navascues, Mohanty & Jayawardhana 2004; Mohanty 
et al. 2004a, b). The remaining 49 are new, comprising 22 sources in Taurus, 5 in IC348, 13 
in Cha I and 9 in Upper Sco. 

Our sample includes most known BDs (i.e., spectral type > M6.5) in the surveyed 
regions, down to Ic ~ 19. In Taurus, Cha I, Upper Sco, TW HyA and R CrA, this includes 
nearly all the sub-stellar objects identified to date (barring a total of ~ 10 objects; of course, 
this is at the time of writing - there are many ongoing surveys for new BDs). In p Oph and 
IC348, our high-resolution optical spectroscopy is limited by a combination of extinction (in 
both) and distance (in IC348). Consequently, we have only 4 sources in p Oph, including 2 of 
the 5 established BDs (Luhman & Rieke 1999); similarly, we could observe only a fourth of 
the ~20 BDs identified so far in IC348 (Luhman et al. 2003a). Moreover, we are somewhat 
biased towards accretors in IC348, in the M5-<M6 range. Luhman et al. (2003a) have found 
~50 objects at these types in IC348; with limited telescope time in hand, we preferentially 
observed those that were likely to be accreting (based on the Ha equivalent widths in Luhman 
et al.'s published low-resolution spectra). As a result, 4 out of our 6 IC348 stars in the M5- 
<M6 range indeed turn out to be accretors; in contrast, the majority of Luhman et al.'s 
objects at these spectral types actually evince low Ha emission (equivalent widths < 10A) 
and are unlikely to be accreting. However, with our IC348 sample contributing a mere 6 
VLMS in this spectral type range to begin with, this bias does not affect our analyses of 
accretor fractions (as we discuss later). Finally, there is no such significant bias in our IC348 
sample at types > M6 (where we have observed two-thirds - 9 out of 14 - of the objects with 
Ic < 19 in Luhman et al.'s study, regardless of whether the latter authors' low- resolution 
spectra suggested possible accretion or not). 

We obtained high-dispersion optical spectra of the sample over numerous observing runs 
from May 2002 to December 2003. For p Oph, Taurus, IC348 and Upper Sco, we used the 
High Resolution Echelle Spectrograph (HIRES; Vogt et al. 1994) on the Keck I telescope at 
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Mauna Kea. With the a slit-width of 1" .15, the 2-pixel binned spectral resolution is R 
33,000. The instrument yielded 15 spectral orders in the 6390-8700 A wavelength region, 
with gaps between the orders, covering a variety of features related to youth and accretion. 
The data were reduced in standard fashion with IDL routines, as described in Basri et al. 
(2000). For the Cha I, R CrA and TW HyA southern targets, we used the Magellan Inamori 
Kyocera Echelle spectrograph (MIKE; Bernstein et al. 2002) on the Baade 6.5m telescope 
at Las Campanas Observatory, Chile. The instrument provided simultaneous spectra over 
~3200-4800A in the blue and ~4800-8800A in the red, with overlapping orders in each 
wavelength range. Most spectra were obtained with a slit- width of 1", giving resolutions of 
R 25,000 and 19,000 in the blue and red respectively. A few spectra in TW HyA were 
obtained with a narrower, 0".7 slit; these observations, and more details on the Magellan 
spectra in general, are given in Mohanty, Jayawardhana & Barrado y Navascues (2003). The 
data reduction was analogous to that employed for the Keck spectra. Finally, we point out 
that all spectra shown in this paper have been smoothed by a 3-pixel boxcar. 



3. Adopted Parameters 

Our analysis of accretion, in subsequent sections, requires estimates for various physical 
quantities, such as effective temperature {T e ff ), mass and line flux. The values we have 
adopted for the pertinent parameters, and the rationale for our choices, are described below. 

Spectral Type: Spectral types for our p Oph, Taurus, IC348 and Chal targets - i.e., for 
most of our sample - are from Luhman & Rieke 1999, Briceno et al. 2002, Luhman et al. 
2003a and 2003b, and Luhman 2004. All these works employ the same optical low-resolution 
SED-fitting technique, with combined dwarf and giant standards as templates, to infer the 
spectral types, ensuring that the latter are all mutually consistent (for p Oph, the types are 
derived in the near-IR, but subsequently pegged to the optical standards; Luhman & Rieke 
1999). For LS-RCrA-1, our only target in R CrA, we adopt a spectral type of M6.5, as found 
by Barrado y Navascues, Mohanty and Jayawardhana (2004) from optical medium- and 
high-resolution spectra. This is consistent with the ~M6 determined for us by K. Luhman, 
with the low-resolution SED-fitting technique (K. Luhman, pvt. comm., 2004). 

For targets in Upper Sco and TW HyA, our adopted types are from Ardila, Martin & 
Basri 2000, Martin, Delfosse & Geieu 2004 (for the USco DENIS targets), and Gizis 2002. 
These authors use various optical molecular band indices for spectral typing, not quite the 
same methodology as SED-fitting. Now, we primarily use spectral types in this paper for 
estimating various physical parameters, such as T e ff and mass, for our accreting objects. For 
a self-consistent analysis, therefore, we must ensure that the spectral types from band indices 
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for our two accretors in Upper Sco and TW HyA - USco DENIS 160603 and 2MASS 1207- 
3932 (§5 and Table 1) - are consistent with the types deduced from SED-fitting in all our 
other accretors. We have checked this by comparing our high-resolution optical spectra of 
these two objects to the spectra of targets with types derived from SED-fitting. Specifically, 
we have examined the TiO bandheads around 8440A , which are highly sensitive to T e ff (Mo- 
hanty et al. 2004a) and should therefore reflect spectral type variations. We find that the 
spectra of USco DENIS 160603 and 2MASS 1207-3932 are consistent with SED-fitting types 
to within ±0.5 subclasses. For instance, 2MASS 1207-3932 has a spectral type of M8 from 
band indices (Gizis 2002); we find that its TiO band shape and strength is comparable to 
that of targets with spectral type M8±0.5, where the latter are derived from SED-fitting. 
Since the intrinsic uncertainty in the SED method is also of order ±0.5 subclasses, we may 
adopt the quoted types for USco 160603 and 2MASS 1207-3932 unchanged. Similar analysis 
of all the non-accreting Upper Sco and TW HyA sources shows that their spectral types, 
from band-indices, agree to within ±0.5-1 subclass with the SED-fitting types. The latter 
variations are somewhat larger than ideally desirable; nevertheless, since we do not use spec- 
tral types for any detailed analysis of the non-accretors, we may adopt their quoted values 
unchanged as well, without affecting the results of this paper 1 . 

Effective Temperature: We estimate T e jj from spectral types via the conversion scale 
specified in Luhman et al. (2003a). This scale is explicitly constructed to agree with the 
predictions of the Lyon evolutionary tracks (Baraffe et al. 1998, Chabrier et al. 2000), which 
are the most widely used tracks today in the low-mass regime. This does not guarantee the 
absolute accuracy of the inferred T e ff : their veracity depends on the validity of the Lyon 
models, which remain somewhat uncertain for very young low-mass objects (Mohanty et al. 
2004a,b). However, the numbers are likely to be in the right ballpark (within 100-200K; 
Mohanty et al. 2004a), as are the relative differences in T e jj over the M5-M9.5 spectral 
subclasses considered here (Mohanty et al., in prep.). Since we are primarily interested in 
the relative variations and trends within our sample, the Luhman-scale is adequate here. 

Age: For the star-forming regions and young clusters in our sample, we adopt ages from 
the literature derived via comparisons to theoretical H-R diagrams. The ages are: <1 Myr 
for p Oph (Luhman & Rieke 1999); ~1.5 Myr for Taurus 2 (Briceho et al. 2002; Luhman 



1 For two sources in our sample - USco 130 and 131 - Ardila, Martin & Basri assign spectral types of 
M7-M8 and M6.5 respectively, from I-J color alone and not low-resolution spectroscopy. Bases on our 
high-resolution TiO analysis, we assign both a spectral type of M7, consistent with their color-based types. 

2 Briceno et al. 2002 state 1-2 Myr; we therefore adopt a mean age of 1.5 Myr. This agrees with the 
<~1.5 Myr age found by White et al. 1999 for the GG Tau quadruple system in Taurus. Luhman et al. 
2003a compare their IC348 sample to that of Brinceho et al.'s low-mass sample in Taurus, and show that 
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et al. 2003a); ~2 Myr for IC348 and Chal (Luhman et al. 2003a, Luhman 2004); ~3 Myr 
for R CrA (Wilking et al. 1997); ~5 Myr for Upper Sco (Preibisch et al. 2002); and ~10 
Myr for TW HyA (Webb et al. 1999; Jayawardhana et al. 1999). For Taurus, IC348, Chal 
and Upper Sco, these mean ages are inferred from the Lyon tracks, while D'Antonella & 
Mazzitelli models (1994, 1997; hereafter DM models) have been used for p Oph, R CrA and 
TW HyA. Ideally, all the ages should be based on the same set of tracks. However, this 
does not affect our statistical analysis of accretion versus age (§5.4), in which we ignore R 
CrA and TW HyA altogether (since we have only 1 accretor in each) . We do include p Oph; 
however, a visual inspection of the Taurus and p Oph HR diagrams (Briceno et al. 2002; 
Luhman & Rieke 1999) suggests that p Oph is indeed somewhat younger than Taurus in an 
absolute sense: sources with similar T e ff are on average more luminous in p Oph. Thus, with 
Taurus at about 1.5 Myr on the Lyon tracks, <1 Myr appears valid for p Oph even with the 
latter models (similarly, though we ignore TW HyA in our age analysis, we note that it is 
surely older than the other regions, regardless of the precise age implied by a set of tracks). 

Mass and Radius: A unique mass can be determined for any source from the theoretical 
evolutionary tracks, given either T e jj and bolometric luminosity (L& /), or T e ff and age. The 
masses of all our IC348 and Taurus accretors (except CHFT-BD-Tau 4) have been supplied 
by K. Luhman (pvt. comm. 2004), and inferred from [T e ff ,-^6oi] comparisons to the Lyon 
theoretical models (the track comparisons are shown in Briceno et al. 2002 and Luhman et 
al. 2003a, though the masses are not explicitly stated there). The T e // used here are from 
spectral types, as described earlier, and L bo i derived by combining J magnitudes, extinction 
estimates and bolometric corrections (Briceno et al. 2002) . For the remaining accretors, we 
have obtained masses by comparing [T e ff , age] to the Lyon models, where 'age' is the mean 
age of the cluster to which the source belongs (given above for the various regions). 

In practice, L hoX or age has a minimal effect on the inferred mass: our cool sources evolve 
down approximately vertical Hayashi tracks at nearly constant T e jj for the first ~10 Myr, 
so T e ff alone is sufficient to derive a relatively precise mass; the L& / or age provides only a 
small correction. This is why we are justified above in using the mean cluster age as a proxy 
for the real age of a given source, even though objects in any specified region scatter around 
this mean. This is also why it is not inconsistent to use Lyon tracks to derive masses for 
accretors in p Oph, Upper Sco and TW HyA, even though the mean ages of these regions 
are from DM models: while it would be technically more consistent to use mean cluster ages 



the Taurus sources appear slightly younger overall than the IC348 ones: they ascribe a median age of ~1 
Myr to Taurus and ~2 Myr to IC348. Our adopted ages for Taurus and IC348 - 1.5 and 2 Myr respectively 
- are commensurate with the statements by Briceno et al., Luhman et al. and White et al., and preserve 
the slight age difference between the two regions found by Luhman et al. 
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also from the Lyon models in these cases, any model-dependent age variations are too small 
to affect the derived masses at these ages once T e ff is specified (and our T e jj in all cases are 
ultimately from the Lyon models, which the Luhman spectral type-T e jj scale is based on). 
As an example of how little mass is affected by the adopted age, we note that for an object 
with T e ff = 2700K (M8 on the Luhman scale), the Lyon models imply a mass of 30 Mj at 
1 Myr and 35 Mj at 10 Myr - a negligible difference for our purposes here. It should be 
pointed out that age does make a difference, in the Lyon tracks, for the very lowest masses 
(< 20 Mj ), which get rapidly cooler and fainter from 1 to 10 Myr. However, we only use 
masses for accretors: among the latter, our two sources at <20 Mj (KPNO-Tau 4 and 12) 
have masses derived via K. Luhman's explicit comparison of [T e ff , L bo [] to the Lyon tracks, 
without making any simplying assumptions about the constancy of T e ff with age. 

The difference between the Lyon and DM models is important when it comes to the 
mass that corresponds to a given T e // . At young ages, and for cool objects like ours, the 
DM models tend to yield masses that are about a factor of 1.5-2 lower than the Lyon ones, 
for a specified temperature (MHCBH03). Insofar as the lowest masses covered in our study 
therefore, our estimates using the Lyon models are more conservative than studies employing 
DM tracks (e.g., MHCBC03, N04). Even with Lyon tracks, however, we go down to ~15 
Mj , i.e., to nearly the planetary-mass (deuterium-burning) limit at 12 Mj . 

Finally with regard to mass, we note that in the Lyon models, the stellar/sub-stellar 
boundary at 0.075 M Q occurs at T e jj m 3000-2930K, over the <1 10 Myr range considered 
here. This corresponds to M6-M6.5 on the Luhman spectral type-T e jj conversion scale, and 
accounts for our conservative working definition of BDs (§1) as objects at types > M6.5. 
While this is accurate for later types, however, one should always keep in mind (as stated 
in §1) that sources close to the boundary could fall either way, due to both spectral typing 
and evolutionary model uncertainties : an object that is M6.5 in our scheme, and whose 
mass is determined via theoretical (Lyon) tracks, can either be a VLMS or a BD, though it 
is guaranteed to lie close to the sub-stellar edge. 

Lastly, accretion rates for some of the higher mass, classical T Tauri stars (CTTs) we 
shall examine have been derived by other groups using veiling measurements, and masses 
and radii from DM models (§5.3). Since M from veiling depends on 7Z*/A4*, we must modify 
these M to Lyon masses and radii, in order to use them in our study. These sources are all in 
Taurus; we have re-derived their masses and radii from Lyon models through the procedures 
outlined above (i.e., converting spectral type to T e ff , and then inferring mass and radius 
by comparing T e jj to the Lyon tracks for the adopted mean Taurus age of 1.5 Myr). Now, 
radius depends more sensitively on age than mass does. However, since our 1.5- Myr estimate 
for Taurus is derived from the same Lyon tracks, not DM ones, our analysis is self-consistent 
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insofar as the theoretical models used. Moreover, the spread in Taurus ages is ~ ±2 Myr; 
over this time and for CTTs masses, the Lyon tracks indicate a variation in (and hence 
in M for a given mass and veiling) of only ~20%. This uncertainty is negligible, since the 
M from veiling are uncertain by factors of ~3-5 to begin with (e.g., Gullbring et al. 1998). 

Equivalent Width: We measure emission equivalent widths for Ha, He I A6678, O I A8446 
and Ca II A8662. Now, in cool mid- to late-M objects such as our targets, no 'true' continuum 
exists in the optical, only a pseudo-continuum comprising a plethora of overlapping molecular 
lines. All our widths are measured relative to the latter, and are thus more precisely pseudo- 
equivalent widths; we drop the 'pseudo-' appellation henceforth only for the sake of concision. 
Furthermore, a perusal of the He I A6678 region in field (Main Sequence, non-accreting) mid- 
to late-M dwarfs reveals a bump at the line position, even in objects with fairly weak activity. 
Whether this is true chromospheric emission in He I , or merely a bump coincidentally 
produced by molecular absorption on either side of the line position, is unclear. We have 
therefore compared the spectra of all our young objects (M5-M9.5) to that of the field dwarf 
Gl 406 (M6), and calculated the He I A6678 equivalent width after excluding the bump at 
this position in the Gl 406 spectrum (this spectral region changes very little, in field dwarfs 
from ~M5-M9, so using a single M6 dwarf template is justified for the level of accuracy 
desired here). While this may underestimate the widths slightly, it precludes our assuming 
He I emission where in fact there is none; we err on the side of caution. Similarly, the 
O I A8446 line sits on top of a TiO absorption band (bandheads at AAA8432, 8442, 8451). 
We have thus calculated O I widths relative to the pseudo-continuum formed by this band, 
through comparison with young, clearly non-accreting objects of the same spectral type in 
our sample with no evidence of O I emission (field dwarfs are not appropriate templates here, 
since the TiO band strengths, though primarily dependent on T e ff , do change in going from 
very young objects at a given spectral type, to field dwarfs of the same type but at much 
higher gravity). Within the context of this methodology, the errors in our inferred equivalent 
widths are about ±10% for widths > lA, and ±0.lA for widths < lA. Our detection limit 
is ~ O.lA. 

Continuum Flux: Our derivation of accretion rates (§5.2, 5.3) involves calculating 
Ca II emission fluxes (J-'caii)- To infer the latter from the observed line equivalent widths 
(EWc a ii), we must know the photospheric continuum flux [T con t) underlying the Ca II lines. 
In previous studies of CTTs, T CO nt has been estimated from the observed /-magnitude. How- 
ever, we shall examine masses all the way from CTTs to BDs, spanning AT e // >1000K. Over 
this range in temperatures, changes in the broad-band /-magnitude reflect a plethora of vari- 
ations in photospheric opacity not necessarily confined to the immediate spectral vicinity of 
Ca II . To circumvent this problem, we use the T CO nt predicted by the latest Lyon synthetic 
spectra for a specified T e // and surface gravity. We use the Lyon STAR-DUSTY 1999 and 
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2000 models 3 (Allard, Hauschildt & Schwenke 2000, hereafter AHS00; Allard et al. 2001), 
T e ff inferred from our spectral types, and a surface gravity of log [g] =4.0 (cgs units). 

The Lyon DUSTY models incorporate the formation of atmospheric dust grains through 
chemical equilibrium calculations, and include both the resulting depletion of grain-forming 
species and dust opacity. These models are appropriate to the mid- to late-M spectral 
types our low-mass sample covers, since field dwarf studies indicate grain formation becomes 
important by late M 4 . They are also appropriate for hotter, higher mass stars, in which no 
dust forms: since dust formation in the models is handled self-consistently, grains simply 
do not appear in the DUSTY models at high enough T e ff . The DUSTY models come in 
two flavors: BD-DUSTY and STAR-DUSTY. Both use the same set of TiO opacities in the 
optical (AMES-TiO line-lists), but different sets of H 2 opacities in the near-IR (AMES-H 2 
line-lists in the BD models and MT-H 2 in the STAR; see AHS00 and references therein). 
Neither flavor is preferred over the other: both BD-DUSTY and STAR-DUSTY have some 
inadequacies in reproducing the observed near-IR SED of M stars, but both match the 
optical colors and SEDs of M stars quite well (AHS00). Under the circumstances, we choose 
the STAR-DUSTY models, for the following reason. In our analysis of accretion rate versus 
Ca II flux (§5.2), we will compare our low-mass, mid- to late-M sample to hotter, higher-mass 
CTTs (mostly K7-M0, but two at K3-K5). Only the STAR-DUSTY models are currently 
available over a large enough range of T e ff (2500-5000K) to simultaneously analyse both 
the high and low-mass regimes. In particular, STAR-DUSTY 1999 models are available up 
to 5000K; STAR-DUSTY 2000 models, which incorporate the latest improvements in dust 
opacity and are therefore better for the mid- to late M's where dust might be important, are 
available up to 4000K. We therefore use the 2000 models for all sources with T e ff < 4000K 
(K7 and later types), and the 1999 ones for the two CTTs with T e// >4000K. 

Our choice of surface gravity is dictated by the mean value indicated by the Lyon 
evolutionary tracks, for low-mass objects that are a few Myrs old. In practice, the continuum 
flux depends predominantly on T e ff , and only negligibly on gravity as long as log[g] ~ 
4.0±0.5 dex (well within the limits predicted by the Lyon models for ages ~1 10 Myrs). 
The T con t adopted is the average value in the synthetic spectra over 8600-8700A , which 
covers the 8662 A Ca II component of interest here (§5). This technique has been profitably 



3 Opcnly available on the internet at ftp://ftp.ens-lyon.fr/pub/users/CRAL/fallard 

4 In practice, dust affects the models, at the wavelengths of the Ca II IRT, only at the lowest T e ff (latest 
types) in our sample: over most of our T e // range, grains have a negligible effect on the Ca II region in the 
synthetic spectra. This agrees with the observation that in field dwarfs, dust appears in blue spectral regions 
at a type of ~M7 and in the red regions appropriate to Ca II at M8-M9, and that in very young objects 
dust appears at even later types than in field dwarfs (see detailed discussion in Mohanty et al. 2004a). 
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used previously for calculating Ha fluxes in field M dwarfs, by Mohanty & Basri (2003); we 
refer the interested reader to that work for further details. Once the continuum flux J- con t is 
known, the line flux is calculated as: Tcaii = ^cont x EWcaii x (1 + r Ca ii)- Here r Ca ii = 
J~ excess I J~cont is the accretion-induced veiling near the Ca II lines, i.e., a measure of the excess 
continuum emission (J- eX cess) produced by the accretion shock relative to the photospheric 
continuum. For objects from other studies, we adopt the quoted rcaii (§5.2); for our own 
sample, we estimate r Ca ii ~ in most cases (§5.3). 

There might certainly be systematic offsets between the fluxes we infer in this manner, 
and the true values, due to uncertainties in the synthetic continuum fluxes. However, this 
effect is likely to be small since, as we will show, our values in the CTT regime agree quite 
well with those derived by other researchers using /-magnitudes. Moreover, we will mainly 
be concerned in this paper with the correlation between Ca II fluxes and accretion rates. 
Once we quantitatively establish this relationship for objects with known M and Ca II fluxes 
found via our technique, we will use it to infer accretion rates for other objects with Ca II flux 
derived in the same fashion. Since the relationship is anchored using M that are calculated 
independent of the Ca II fluxes, the "true" value of the flux is not key to the analysis; it is 
only important that all the fluxes be derived via a consistent methodology. 

Rotation Velocity: We derive rotation velocities (v sini ) of our targets by cross-correlating 
with a "spun-up" template of a slowly rotating standard. Multiple spectral orders (~ 6), 
selected on the basis of an absence of strong telluric features, strong gravity-sensitive features 
and stellar emission lines, were used for the analysis. Following White & Basri (2003), we 
used a combination of giant and dwarf spectra for our templates. A detailed discussion of 
the cross-correlation technique is given in Mohanty & Basri (2003). Our errors are, conser- 
vatively, ~ ±2 km s _1 for v sim < 20 km s -1 , ±2.5-3 km s -1 for v sim of 20-50 km s _1 , and 
±5 km s _1 for v sim of 50-70 km s _1 (~60km s _1 being the highest v sim in our sample). 



4. Accretion Diagnostics 

It is useful to begin with a summary of the expected behaviour, in the presence of 
accretion, of the main diagnostic spectral lines included in our spectra: Ha (6563A), He I 
(6678A), O I (8446A) and Ca II (8662A, the reddest component of the Ca II infrared triplet). 



4.1. Ha 



In classical T Tauri stars (CTTs), the widely accepted paradigm of magneto-spheric 
accretion holds that stellar field lines threading the disk channel material from the inner 
edge of the disk onto the stellar surface. The gas in these 'accretion columns' falls in at 
nearly free-fall velocities, resulting in an accretion-shock as it strikes the surface; Ha is 
produced both in the shock region and in the infalling gas. Consequently, Ha emission is 
generally much stronger in CTTs than in non-accreting weak-line T Tauri stars (WTTs), 
where Ha arises only from chromospheric activity. The Ha line-profiles in CTTs are also 
distinctive, evincing broad line-wings extending out to a few hundred kilometers per second 
from line-center; the broadening is due to both the high (free-fall) velocity of the accreting 
gas and the large damping wings associated with Stark-broadening in this very optically 
thick line (Muzerolle, Calvet & Hartmann, 2001). Often, the CTTs Ha emission is also 
asymmetric, due to inclination effects (e.g., the star/disk may occlude part of the infall) 
and/or absorption by an accompanying outflowing wind (e.g., Alencar & Basri 2000). 

The strength of the Ha line has long been used as a simple basis for distinguishing 
between CTTs and WTTs, with objects exhibiting Ha equivalent widths > 10A classified 
as CTTs. However, this diagnostic is not applicable over any large range of spectral types, 
due to the 'contrast effect'. Basically, as one moves to cooler (later type, lower mass) objects, 
the same Ha flux stands out better against the increasingly faint photosphere. As a result, 
the equivalent width of even chromospheric Ha can exceed 10A at late spectral types. 
As an alternative, White & Basri (2003; hereafter WB03) have proposed a quantitative 
discriminant between CTTs and WTTs based on Ha line profiles. In CTTs, the accretion 
shock produces excess continuum emission in the optical that fills in underlying photospheric 
absorption lines; such 'veiling' is generally a sure sign of accretion. WB03 point out that, 
empirically, the full-width of the Ha line at 10% of the line-peak invariably exceeds 270 km 
s _1 in veiled CTTs, regardless of spectral type. Thus an Ha 10% width > 270 km s _1 can 
by itself be adopted as a good accretion indicator, without a laborious check for veiling. 

Subsequently, JMB03 concluded that WB03's accretion criterion, based primarily on 
CTTs observations and linked explicitly to the presence of veiling, is too restrictive when 
applied to the VLMS and BD regimes. Accretion rates in the latter (~ 10~ 10 -10~ 12 M 
yr _1 ) are generally 1-3 orders of magnitude lower than in CTTs, too low for any appreciable 
veiling (Muzerolle et al. 2000, 2003). At such puny rates, Stark-broadening of the Ha line 
is also negligible; the line-broadening is dominated by the infall velocities alone, and thus 
smaller than in CTTs. Furthermore, VLMS and BDs are less massive than CTTs, implying 
lower free-fall velocities for the gas and thus a further reduction in the extent of the Ha 
line-wings. From an examination of young low-mass objects in IC348, JMB03 showed that 
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other accretion signatures, such as strong emission in the permitted lines of He I , O I and 
Ca II and asymmetries in the Ha line profile, are associated with Ha 10% widths as low as 
200 km s _1 . They proposed, therefore, that an Ha 10% width > 200 km s _1 is a better 
accretion diagnostic for VLMS and BDs. 

JMB03's suggested accretion cutoff is also supported by recent infrared spectroscopy. 
N04 show that young VLMS and BDs with Ha 10% width > 200 km s _1 generally also 
exhibit Pa/9 emission associated with accretion. Moreover, they find that in objects above 
this cutoff, the Ha 10% width correlates very well with the mass accretion rate derived 
through other means. In the present paper, therefore, we adopt the 200 km s _1 Ha criterion 
as our main accretion diagnostic (subject to the caveats discussed below). This will be 
supported by line-profile asymmetries, as well as emission in other accretion-related lines 
(§4.2) where applicable. 

Given that the infall rates in young VLMS and BDs are generally very small, and that 
the Ha line-wings are likely to be dominated by infall velocities, three caveats must be borne 
in mind when using the Ha line profile as an accretion diagnostic. First, the rotation velocity 
of the object becomes important. As we will show, the intrinsic full-width of the Ha line in 
chromospherically active, non-accreting young VLMS and BDs in our sample is usually ~ 
100-150 km s _1 . In a VLMS or BD with v sini >50 km s" 1 , the full-width of the Ha profile 
will be rotationally broadened by a further >100 km s _1 . Thus, with an adopted accretion 
cutoff of 200 km s" 1 , a merely active but rapidly rotating low-mass object can appear to be 
accreting. This effect is indeed present in a few of our targets, as we will discuss in §5.1. 

Second, since it is not the absolute infall velocities but rather the apparent ones along 
the line of sight that produce the Ha line-wings, inclination effects are important - more 
so than in CTTs, where inclination-independent Stark-broadening usually dominates in the 
wings. In an accretor oriented pole-on to our line of sight, most of the high- velocity infalling 
gas moves along the plane of the sky; hence the emission is concentrated at line-center (e.g., 
Muzerolle et al. 2003; hereafter MHCBH03). Thus, in a weakly accreting VLMS or BD seen 
pole-on, the Ha profile will be narrow and symmetric, possibly with a 10% width < 200 km 
s _1 , reminiscent of chromospheric activity instead of accretion. This effect too is apparent 
in our sample, and must be corrected for by examining other accretion-related lines (§5.1). 

Finally, the observed Ha 10% width depends on the line optical depth. At very low 
accretion rates, there simply may not be enough infalling gas to push the width above our 200 
km s^ 1 cutoff. An examination of the Ha models profiles and comparisons to observations 
shown in MHCBH03 indicates that this effect is likely to be important for mass infall rates 
< 10~ 12 M Q yr _1 ; of course, it is exacerbated by the inclination-angle effect discussed above. 
In summary, adopting an Ha 10% width > 200 km s _1 as an accretion indicator is probably 
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not sufficient to identify the weakest accretors (M < 10~ 12 M yr -1 ). However, it is good for 
a conservative estimate of the number of accretors in our sample, especially when combined 
with the other accretion diagnostics addressed in the next section. 

Here we would like to clarify the physical basis for JMB03's 200 km s _1 accretion cutoff 
for VLMS and BDs. For masses < O.1M , free-fall velocities from infinity are of order 200 
km s^ 1 , which would simplistically lead to Ha full-widths of ~400 km s _1 . However, the 
gas is expected to flow in from the inner edge of the disk; a variety of studies indicate that 
this lies at ~ 2-3 sub-stellar radii in the VLMS and BD regime (e.g., Mohanty et al. 2004c). 
Free-fall velocities from these radii are only ~150 km s" 1 ; since the most likely inclination 
angle for a random distribution of orientations is ~60°, the apparent infall velocity is further 
reduced on average to ~125 km s" 1 , implying an Ha full-width of ~250 km s _1 . Moreover, 
since we consider the full-width at 10% of the peak, and not at the true 'base' of the line, we 
measure less than the full range of infall velocities (e.g., the width at say 5% is larger; 10% is 
an ad hoc value chosen to ensure a good measurement of the width without confusion by the 
photospheric continuum). Finally, since the line- width depends on the optical depth, which 
decreases as one moves to larger offsets from the mean velocity, we will always measure a 
lower limit to the true range of gas velocities (an effect that is exacerbated with declining 
accretion rate). These considerations indicate that a lower limit of 200 km s _1 in Ha 10% 
width is a physically reasonable accretion criterion. As JMB03, MHCBH03 and N04 have 
shown, other accretion indicators are indeed present in objects with Ha 10% widths larger 
than 200 km s -1 , so this limit is empirically justified as well. 



4.2. He I , O I and Ca II 

Permitted emission lines of He I , O I and Ca II are also associated with accretion 
in CTTs (Muzerolle, Hartmann & Calvet 1998, hereafter MHC98) - in particular, He I 
AA5876 and 6678, O I AA7773 and 8446, and the Ca II infrared triplet (IRT) lines AAA8498, 
8542 and 8662. Of these, our spectra include He I A6678, O I A8446 and Ca II A8662. In 
their detailed study of CTTs, MHC98 showed that: the O I emission, as well as the broad 
(FWHM >100 km s _1 ) components of the Ca II IRT and He I lines, probably originate in 
the magnetospheric infall region; O I A8446 appears in emission only in CTTs with accretion 
rates (M ) > 10~ 8 M yr -1 ; and He I , O I and Ca II line strengths all increase with M 
in CTTs, with the Ca II IRT exhibiting the most striking correlation. In young VLMS and 
BDs, the recent analyses by JMB03 and MHCMH03 reveal the following. 

I : JMB03 find O I A8446 in some, though not all, of their low-mass IC348 accretors 
(defined as objects with Ha 10% width > 200 km s _1 ); the emission profile is broad, with 
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FWHM > 100 km s" 1 . None of their non-accretors shows this line in emission. MHCBH03 
do not report on the behavior of permitted O I . However, by modeling the Ha line profile, 
they have derived infall rates for objects in their sample, including some of the accretors in 
which JMB03 detect O I A8446. The inferred M are of order 10~ 10 M yr _1 , showing that 
in VLMS and BDs, this line can appear in emission at much lower accretion rates than in 
CTTs. This is not surprising: a weak emission feature, that would be invisible against the 
bright photospheres of CTTs, is expected to show up better against the intrinsically faint 
photospheres of much cooler low-mass objects (the 'contrast effect' mentioned earlier). 

He I : MHCBH03 find that emission in He I A5876 is ubiquitous, in both accreting 
as well as chromospherically active VLMS and BDs. They do not discuss the behavior of 
He I A6678 in their sample. JMB03 do find He I A6678 emission in one IC348 accretor 
(IC348-165), and do not report any emission in their non-accretors 5 . 

Ca II : MHCBH03 find widespread emission in the Ca II IRT as well. A closer perusal 
of their results shows that the AA8498 and 8542 IRT components are present in both accret- 
ing and merely active objects; however, the A8662 component occurs almost exclusively in 
accreting objects (9 of their 11 A8662 detections come from accretors, though not all their 
accretors exhibit this line). Similarly, JMB03 detect Ca II A8662 in some but not all of their 
IC348 accretors; in two of their three detections, the emission is relatively broad (> 100 
km s _1 ) 6 . Crucially, and in agreement with MHCBH03, they do not find any Ca II A8662 
emission in the non-accretors. 

Thus, emission in O I A8446, He I A6678 and Ca II A8662 seems reasonably well- 
correlated with accretion in VLMS and BDs, as in CTTs: while not detected in all the low- 
mass accretors, these lines appear preferentially associated with accretion, and not activity. 
This assertion seems particularly safe to make for Ca II A8662, for which the combined size 
of the JMB03 and MHCBH03 samples is large. We conclude that emission in He I A6678, 
O I A8446 and Ca II A8662 is suggestive of accretion, and certifies accretion when allied 



5 A review of the spectra we presented in that paper shows He I A6678 emission is also present in IC348- 
265, as noted in Table 1. With an Ha 10% width slightly less than 200 km s _1 , this was not classified as 
an accretor by JMB03; however, we will argue in §5 that this object is in fact likely to be accreting as well. 

6 MHCBH03 state that broad emission in O I , He I and Ca II is present only in the earliest type 
accretor in their sample, LkHa 358 (K7-M0), and not in any of their later- type accreting objects. This is 
a little surprising at first, since JMB03 point out broad emission in some of the same low- mass objects in 
MHCBH03's sample. The problem lies in the definition of "broad": MHCBH03 define this as a FWHM > 
200 km s _1 , while JMB03, and we in this paper, define it as a FWHM > 100 km s _1 , in order to remain 
consistent with the definition in MHC98. With this defining value, it is likely that O I , He I and Ca II 
emission from some of MHCBH03's accretors is broad as well, at least in targets in common with JMB03's. 
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with an Ha 10% width > 200 km s 1 . We will re-examine this issue with our data in §5. 
In contrast, based on MHCBH03's results, He I A5876 and Ca II AA8498 and 8542 are not 
good accretion indicators in the low-mass regime, appearing in emission in both accreting 
and active objects. We will not discuss the latter lines any further. 

5. Results and Discussion 

5.1. Observed Accretion and Ouflow 

Our accretion and activity results are summarized in Table 1. Ha profiles for half the 
sample have been shown in our previous work (JMB02; JMB03; Mohanty, Jayawardhana & 
Barrado y Navascues 2003; Barrado y Navascues, Mohanty & Jayawardhana 2004); profiles 
for the new objects in this paper are plotted in Figs. 1 and 2 (all accretors, old and new, 
in Fig. 1; new active objects in Fig. 2). The majority of our 82 targets evince narrow and 
symmetric Ha lines, with 10% widths ~ 100-150 km s _1 (Fig. 2). Such emission is expected 
from chromospheric activity alone, so we classify these objects as non-accretors. 19 objects 
in our sample, however, have Ha 10% widths > 200 km s" 1 (Fig. 1), as well as low rotation 
velocities (v sim < 20 km s _1 ). Given the small implied rotational broadening, and the much 
narrower profiles in the sources that are merely active, the very broad emission in these 19 
targets is highly unlikely to result from a combination of activity and rotation. Instead, such 
profiles are expected for gas accreting at nearly free-fall velocities (§4.1); we thus classify 
these sources as "bona-fide" accretors. Asymmetries in many of their Ha profiles, analogous 
to those often observed in CTTs, support this conclusion. 

A further two objects, IC348 355 (spectral type M8) and USco 75 (M6), also have 
Ha 10% widths slightly above 200 km s -1 ; however, they are very rapid rotators, with 
a v sim of 45 and 63 km s _1 respectively. Their large Ha 10% widths are thus likely to 
result primarily from rotational broadening, not accretion. This conclusion is strengthened 
by their Ha equivalent widths (which are unaffected by v sim effects): these are relatively 
low, and comparable to those in other non-accretors of the same spectral type. They also 
do not evince any emission in He I , O I or Ca II , unlike many of the accretors (see further 
below). We thus classify these two objects as non-accretors 7 . 

Finally, we note that while MHO-5 shows no signs of accretion in our spectrum - its Ha 
10% width is less than 200 km s" 1 , and it has no emission in permitted He I , O I or Ca II 
- MHCBH03 have detected strong forbidden [O I ] emission at AA6300 and 6363 from this 



7 Note that JMB03 classified IC348-355 as an accretor; however, they did not consider rotation effects. 
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object. Since such emission is associated with mass outflows in CTTs, which in turn provide 
indirect evidence of ongoing accretion, MHCBH03 classify MHO-5 as an accretor, though 
its Ha 10% width is less than 200 km s _1 in their data too. Our Keck spectrum of MHO-5 
does not cover these forbidden O I lines, so we cannot confirm MHCBH03's result; under the 
circumstances, we simply adopt their evaluation of this object as a bona-fide accretor. In 
this case, its relatively narrow Ha profile is probably due to a pole-on viewing orientation, 
as MHCBH03 also point out. 

Now let us examine the prevalence of He I , O I and Ca II emission. Out of the 15 
targets in which He I A6678 is detected, 11 are likely accretors by our Ha analysis above. 
Similarly, 13 of the 17 Ca II A8662 detections are in accretors, as are all 12 of the O I A8446 
detections (Ca II emission shown in Fig. 3). This strongly validates our previous conclusion, 
based on the results of JMB03 and MHCBH03, that emission in these lines correlates with 
ongoing accretion. In turn, this result prompts us to re-examine the status of the few objects 
which have Ha 10% widths < 200 km s _1 , but nevertheless do evince emission in at least 
one of these lines: GY 310, KPNO-Tau 4, KPNO-Tau 11, CFHT-BD-Tau 4, MHO-4, IC 348 
256 and USco 128. 

Table 1 shows that GY 310 and IC 348 256 in fact have Ha 10% widths reasonably 
close to our accretion cutoff of 200 km s _1 - 144 and 180 km s _1 respectively, comparable to 
the 166 km s _1 in the apparent accretor MHO-5 discussed above. Fig. 1 moreover reveals 
that the Ha profiles of the two sources are somewhat asymmetric. Finally, there is direct 
evidence for inner disks in both, from mid-IR excesses (Mohanty et al. 2004; Jayawardhana 
et al. 2003). While the presence of a disk does not certify accretion, since circumstellar 
material may exist beyond the main accretion phase (§5.5), it is certainly consistent with 
ongoing disk-accretion. These factors, combined with emission in either Ca II or He I , leads 
us to classify GY310 and IC348 256 as "probable" accretors. 

Similarly, the Ha 10% widths in KPNO-Tau 4 and 11 and CFHT-BD-Tau 4 are not 
much less than 200 km s" 1 (155, 175 and 182 km s -1 respectively). However, their Ha profiles 
appear quite symmetric (Fig. 1), similar to those in chromospherically active objects. We 
therefore classify them as "possible" accretors, with the following caveats. We have already 
seen that Ca II A8662 emission is associated only with accreting objects in the rest of 
our sample, implying that any purely chromospheric component of this line is very small. 
Nevertheless, even such a weak component may stand out in emission against the faint ultra- 
cool photosphere of KPNO-Tau 4 ('contrast effect'); at M9.5, this is the coolest object in 
our sample. Consequently, a purely chromospheric origin for its Ca II emission is not wholly 
implausible. On the other hand, Mohanty & Basri (2003) have shown that such activity falls 
off sharply, at least in field dwarfs, at ~M9, with little or no Ha emission at later types; 
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Mohanty et al. (2002) have further argued that this fall-off is due to the increasing neutrality 
of the cool atmospheres with later type. Whether the same applies to young objects, i.e., 
whether a very cool young BD like KPNO-Tau 4 can in fact sustain sufficient activity to 
produce the observed Ha and Ca II emission, is unclear. Spectroscopy of a large sample of 
young ultra-cool objects (very few known so far) is needed to settle this question. 

In CFHT-BD-Tau 4, which at M7 is somewhat hotter than KPNO-Tau 4, the contrast 
effect is not as strong. Nevertheless, it is still quite cool, and its Ca II emission is also 
very weak, so chromospheric Ca II cannot be completely ruled out. It is known to have a 
disk from mid-IR and sub-mm/mm observations; however, the disk also appears relatively 
evolved (Sterzik et al. 2004). These facts do not permit a clear distinction between activity 
and accretion, so 'possibly accreting' is the most conservative conclusion. Finally, in KPNO- 
Tau 11, the contrast effect is much less important in both above sources, since the latter is 
siginficantly hotter at M5.5. However, KPNO-Tau 11 does not show any Ca II A8662, only 
He I A6678. Though this line is also predominantly associated with accretion in our sample, 
we do find it, without any accompanying Ca II A8662 emission, in two objects we classify 
as non-accretors (see below). Thus KPNO-Tau 11 may not be accreting, only much more 
active than most of our sample. Further observations are required to clarify its status. 

Lastly, MHO-4 and USco 128, while exhibiting clear He I A6678 emission, evince Ha 10% 
widths significantly lower than 200 km s _1 (115 and 111 km s _1 ), as well as highly symmet- 
ric Ha profiles (Fig. 1). We thus classify them as non-accretors; note that MHCBH03 have 
independently found these two sources to be non-accretors as well. Presumably, their He I 
A6678 emission arises from very strong chromospheric activity, though we caution that accre- 
tion cannot be completely ruled out: very low infall rates combined with severe inclination 
effects may conceivably produce the narrow and symmetric Ha profiles we observe. 

Two further sources deserve mention: KPNO-Tau 14 in Taurus, and DENIS 161929 in 
Upper Sco. While deriving v sini (§3), we discovered that their cross-correlation functions are 
slightly double-peaked, a telltale signature of spectroscopic binaries (SB2s). We are carrying 
out follow-up observations to confirm this finding; for the moment, we refrain from citing 
a v sim for them, and classify them as possible SB2s. Both have Ha 10% widths less than 
200km s" 1 , and no emission in He I , O I or Ca II . A close look at their Ha profiles does 
reveal some asymmetry; e.g., the profile of DENIS 161929 (Fig. 2) is similar to that of the 
"probable" accretors GY310 and IC348 256 (Fig. 1). However, such asymmetries can also 
result from overlapping Ha profiles in an SB2. Without more definitive accretion signatures 
(broad Ha, or emission in other accretion-related lines as in GY310 and 1C348 256), we 
conservatively classify them as non-accretors for now. As an aside, we note that KPNO-Tau 
14 has a spectral type of M6, and DENIS 161929 of M8; if their SB2 status is confirmed, 



-20- 



they would be the two latest-type (lowest mass) spectroscopic binaries discovered to date at 
ages of a few Myrs. As such, they might prove very useful for testing theoretical evolutionary 
tracks for very young VLMS and BDs, through empirical mass and radius determinations. 

In summary a total of 25 VLMS and BDs in our sample evince accretion signatures 
of some sort or another (broad Ha and/or emission in other accretion- related lines). Of 
these, 20 are "bona-fide" accretors: 19 based on Ha profiles, and 1 (MHO-5) based on its 
forbidden [O I ] emission in MHCBH03's data. The remaining 5 sources (GY 310, IC348 
256, KPNO-Tau 4, KPNO-Tau 11 and CFHT-BD-Tau 4) are classified as "probable" or 
"possible" accretors; further observations are required to verify their accretor status. 

Recently, Barrado y Navascues & Martin (2003; hereafter BM03) have suggested equiv- 
alent width criteria for discriminating between infall and activity in both higher-mass stars 
as well as in VLMS and BDs. They set the upper limit for chromospheric Ha flux at 
\og[Lual 'Lboi] = -3.3, the average saturation value observed in young open clusters. This 
flux limit is then converted into a quantitative equivalent width cutoff for different spec- 
tral types. Objects with Ha widths above the cutoff appropriate to their spectral type are 
assumed to be accreting, while those below the limit are likely to be merely active. Note 
that since the same Ha flux becomes more prominent against a cooler photosphere (contrast 
effect), the limiting equivalent width derived by BM03 increases with later spectral type. 

The BM03 scheme is specifically constructed for low-resolution spectra, which generally 
yield larger equivalent widths than high-resolution ones. Later in the paper (§5.4), we will use 
previously obtained low-resolution data for our low-mass sample to derive accretor fractions 
based on the BM03 criteria. It is useful to briefly discuss here, however, the results we will 
obtain. In general, we find that the BM03 scheme performs reasonably well: targets that we 
find to be accreting are largely classified as such by BM03's technique, and similarly for non- 
accretors. Nevertheless, the BM03 scheme systematically gives larger accretor fractions for 
our low-mass sample, when applied to low-resolution Ha widths, compared to the fractions 
suggested by our detailed high-resolution analysis above. There are two reasons for this. 
First, sources that are more active than the norm for their spectral type, or are possibly 
chromospherically flaring, will be mis-classified as accretors by the BM03 criteria; their 
non-accretor status can only be deciphered from a high-resolution examination of their Ha 
profiles such as undertaken here. Conversely, weak accretors seen pole-on can have narrow 
and symmetric Ha profiles at high-resolution, while their Ha equivalent widths and fluxes 
(unaffected by inclination) can still be higher than that in active objects. Such objects would 
appear to be merely active from our high-resolution analysis, while they would correctly be 
classified as accretors by BM03. MHO-5 is a good example of this. Its Ha equivalent width 
at low resolution is substantially higher than in chromospherically active sources at the same 
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spectral type (M6), implying ongoing accretion by the BM03 criteria. However, the object 
appears a non-accretor in our high-resolution data - its Ha profile is narrow and symmetric, 
and it lacks other accretion-related lines such as Ca II . We are aware of its true accretor 
status only serendipitously, from the [O I ] emission spotted by MHCBH03; other MHO-5 
analogs may well lack such emission, and thus simply appear active in our high-resolution 
study. In summary, our high-resolution analysis is more conservative in identifying accretion; 
while we may classify some accretors as active, we are unlikely to deduce accretion in merely 
active sources. While the BM03 scheme is very useful for statistical analyses of large samples, 
we prefer our conservative approach for the classification of individual sources. 

Accretion rates are generally much lower in the low-mass regime than in CTTs (see §5.2 
and 5.3 below); much of the preceding discussion underscores the difficulties that occasionally 
arise in distinguishing between infall and activity under these circumstances. Nevertheless, 
out of the 25 sources in which we have identified some accretion signatures, 20 can be clearly 
classified as bona-fide accretors. The latter sample includes: (1) an accretor at every spectral 
subclass from M5 to M9, (2) the lowest mass accreting BD known to date - the Taurus M9 
object KPNO-Tau 12 (mass ~ 20 Jupiters), (3) the only known BD accretor in the 5 Myr- 
old Upper Sco region, where most stars have also stopped accreting - DENIS 160603 (M7.5, 
mass ~ 40 Jupiters), and (4) the oldest accreting BD currently known - the 10 Myr-old 
TW HyA member 2MASS 1207-3932 (M8, mass ~ 35 Jupiters). Our observations of this 
last object have been discussed in detail in Mohanty, Jayawardhana & Barrado y Navascues 
2003. Our diagnosis of accretion in that paper, based on Ha and He I emission, has since 
been buttressed by X-ray analyses (Gizis 2004), and the detection of accretion shock-induced 
UV emission (Gizis 2004, pvt. comm.). We note that this BD is now also known to possess 
a surrounding disk, from mid-IR excess measurements (Sterzik et al. 2004). 

Next, we turn our attention to mass outflows, which are often associated with accretion 
in CCTS, and are widely identified in the latter through forbidden line emission. We have 
already mentioned that MHO-5 shows strong [O I ] emission in MHCBH03's data, though 
this region is not covered in our Keck spectrum. We do not see any other optical forbidden 
lines (e.g., [N II ], [S II ]) in this object. Similarly, Natta et al. 2003 see hints of [N II ] and 
[S II ] emission in Cha Ha 2 and 6 (which we classify as accretors based on Ha), and in Cha 
Ha 3 (classified here as a non-accretor); we do not see these lines in our spectra of these 
three objects (possibly due to lower S/N). However, we do find very strong emission in [O I ], 
[N II ] and [S II ] in our Magellan spectra of the accretor LS-RCrA-1, as discussed in Barrado 
y Navascues, Mohanty & Jayawardhana 2004 (see also Fernandez & Comeron 2001). At a 
spectral type of ~M6.5, this object resides at the stellar /sub-stellar boundary, and to date 
has offered the best evidence for the existence of outflowing jets and winds in the low-mass 
regime. We have now detected, for the first time, weak [O I ] emission in 2MASS 1207-3932 



-22 - 



as well (Fig. 4). At an age of 10 Myr, and a spectral type of M8 (implying a mass ~ 35 
Jupiters), this TW HyA member is the oldest, lowest mass BD known to evince any mass 
outflow signature (as well as the oldest known BD with accretion signatures, as discussed 
above). The overall paucity of jet/wind detections in VLMS and BDs so far, compared to 
the CTT regime, is not surprising: outflow rates are expected to scale with infall rates, which 
are very small in the low-mass domain. Nonetheless, our detections, combined with those 
by others, imply that such outflows can exist in accreting VLMS and BDs. 

Taken together, the above results suggest that: (1) CTT-like accretion extends over 
the entire sub-stellar regime, down to the lowest mass BDs known, (2) outflows can also be 
present in accreting VLMS and BDs, as in CTTs, and (3) the accretion timescales in VLMS 
and BDs are comparable to those in higher-mass stars (as discussed in detail in §5.4). 

5.2. Accretion Rates 

In CTTs, accretion rates are usually determined through veiling measurements. Unfor- 
tunately, infall rates in VLMS and BDs are usually far too small to produce any discernible 
veiling (WB03, MHCBH03). For these objects, M may be estimated through detailed mod- 
eling of the Ha line profile (MHCBH03). This process is laborious and time-consuming, 
however, and a more straightforward technique is highly desirable. A possible solution is 
to employ Ca II A8662 fluxes: as discussed earlier, Ca II IRT component fluxes correlate 
very well with M at least in higher- mass CTTs (MHC98). In our low-mass sample, we have 
already shown that the presence of emission in the A8662 component is closely allied to on- 
going accretion. We now demonstrate that the flux in this line also correlates quantitatively 
with M for these low masses, just as in CTTs. For low-mass accretors with Ca II emission, 
this will permit a direct derivation of M without resorting to Ha profile-modeling. 

MHC98, WB03 and MHCBH03 provide Ca II A8662 equivalent widths, veiling and M 
measurements for a number of objects, from CTTs to BDs. The accretion rates are derived 
from either the observed veiling (all the CTTs) or detailed Ha line-profile modeling (most 
of the VLMS and BDs). We use this sample to probe the relationship between Ca II flux 
and M from the stellar to sub-stellar regimes. Our analysis technique is as follows. 

(1) We convert the Ca II equivalent widths quoted by MHC98 and MHCBH03 to fluxes 
through the technique described in §3. In the case of CIDA-1 and GM Tau, WB03 provide 
M and veiling but not the Ca II equivalent widths; however, we have access to the spectra 
they used, and derived the Ca II widths in these ourselves. (2) The veiling M we use are from 
MHC98 and WB03; these are all calculated using the methodology outlined by Gullbring et 
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al. 1998. This is explicitly true of the WB03 values (derived from Rc-hand veiling), and 
of most of the MHC98 values, taken directly from Gullbring et al.'s work (who use [/-band 
veiling). For a few heavily- veiled stars not analysed by Gullbring et al., MHC98 use M from 
Hartigan, Edwards & Ghandour 1995 (calculated from V-band veiling), whose assumptions 
leads to systematically higher M than the Gullbring et al. technique; however, MHC98 have 
adjusted these to the Gullbring et al. scale as well, through multiplication by an appropriate 
factor. (3) All these M inferred from veiling are proportional to the assumed 7Z*/M.*, so we 
must ensure that the masses and radii used are consistent with the Lyon tracks adopted in 
this paper. WB03 use Lyon tracks, so we adopt their M unchanged; the M in MHC98 are all 
based on DM tracks, so we have modified these M to Lyon masses and radii (methodology 
given in §3). (4) The M from Ha profile-modeling, on the other hand, are quite insensitive 
to mass and radius assumptions; thus, we adopt unchanged the M found via such modeling 
by MHCBH03. (5) We have left out unresolved binaries included in the MHC98 and WB03 
studies, since binarity can vitiate the calculation of M from veiling (Gullbring et al. 1998). 

The sources we examine, and their final adopted parameters, are listed in Table 2. The 
resulting relationship between Ca II A8662 flux and accretion rate is shown in Fig. 5a. The 
two quantities clearly correlate remarkably well, over more than 4 dex in both M and line 
flux. MHC98 found such a relationship in CTTs, which we recover in Fig. 5a for M between 
~1CT 6 and 1CT 8 - 5 M yr _1 . In particular, Fig. 9 in MHC98 shows the correlation between 
their M and Ca II A8542 fluxes, where the line equivalent widths are converted to fluxes 
using /-band magnitudes and veiling estimates. We have used their veiling measurements, 
as well as their quoted A8662 widths, which are very similar to their A8542 ones; thus, any 
differences between their Fig. 9 and our Fig. 5a for the same objects are ascribable to our 
use of synthetic continuum fluxes to infer line fluxes, versus their use of /-magnitudes for 
the same purpose. A comparison of the plots shows that our results are in fact very similar; 
while our derived Ca II fluxes are systematically slightly higher than theirs, the trend in flux 
versus M is identical (note that a few outliers in their plot are not included by us, since 
they are binaries). This bolsters our confidence in the technique we have used to derive the 
line fluxes. Compared to MHC98, our crucial new result is that the correlation extends all 
the way to VLMS and BDs, with accretion rates several orders of magnitude lower than in 
CTTs (most objects plotted with M < 10" 8 - 5 M yr* 1 are VLMS or BDs; see Table 2). A 
linear fit to the data yields: 

log(M) = 1.06 log(F C aii) ~ 15.40 CTTs + VLMS + BDs [1] 

While a correlation between M and Ca II flux is obvious, however, Fig. 5a also reveals 
significant scatter around our derived fit: the la scatter in M around the fit is a factor of 
3, and the largest deviations are up to factor of 6. Now, emission lines and veiling in CTTs 
are known to vary over time, and variability is also observed in low-mass accretors (see 
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discussion in N04). However, it is unlikely that the observed scatter is due to such variablity 
combined with a difference in epoch between the Ca II flux and M measurements. For all 
the low- mass objects (GM Tau, CIDA-I and all IC348 sources), the Ca II fluxes and M have 
been derived from the same spectra. Similarly, while spectra from different epochs have 
been used for the higher-mass CTTs (Ca II measurements by MHC98, M from Gullbring et 
al. 1998 and Hartigan, Edwards & Ghandour 1995), MHC98 demonstrate that there is no 
significant change in the sources over this time (veilings found by MHC98 correlate very well 
with the veiling-dependent M derived in the earlier studies). In fact, a closer perusal of Fig. 
5a shows that the scatter is systematic: all the low-mass objects lie somewhat above our fit, 
while most of the higher-mass CTTs fall somewhat below. We illustrate this in Fig. 5b, by 
obtaining separate fits to the high- and low-mass regimes; the linear fits derived are: 

log (M) = 0.71 log (F C aii) - 12.66 CTTs [2] 

log (M) = 0.93 log (Fcaii) ~ 15.03 VLMS + BDs [3] 

We immediately see that the scatter around the individual fits is now much reduced: the la 
scatter in M is a factor of ~1.5 for each fit, and the largest deviations are at most a factor 
of ~3. The fits suggest an offset in the Ca II flux-M relationship between the high and 
low-mass regimes, with the low-mass objects evincing more Ca II flux at a given M . 

The reasons for this displacement, however, are unclear. Systematic problems in the 
synthetic continuum fluxes, in going from the hotter CTTs to the cool low-mass sample, do 
not see responsible: the same offset appears even if we use observed J-band magnitudes to 
calculate the Ca II fluxes of the high- and low-mass samples (not plotted). Overestimation of 
the Ca II flux in the low-mass accretors, through contamination by chromospheric Ca II , is 
also unlikely: as discussed earlier, clearly non-accreting low-mass objects in our study simply 
do not show any Ca II emission, so the purely chromospheric contribution is likely to be 
negligible in the low-mass accretors. One might postulate, instead, that there are differences 
between the high- and low-mass regimes in the surface covering fraction of the accretion hot 
shock: if the shock covers relatively more surface area in the low-mass sources for a given 
M , then their accretion-related Ca II flux might also be higher than in CTTs, as observed. 
However, this is not supported by the data: using the observed veilings, estimating the 
underlying photospheric fluxes (e.g., from synthetic spectra), and assuming a shock T e // ~ 
10 4 K (e.g., MHCBH03), one easily finds that the covering fraction in the low-mass sources 
is <Cl%, which is in fact quite similar to that in CTTs (as also shown by MHCBH03). 

One possibility is a mismatch between the M from the ultraviolet [/-band veiling (most 
of the higher-mass CTTs) and M from optical diagnostics (all the low-mass objects: M from 
optical veiling in GM Tau and CIDA-1, and from modeling optical Ha profiles in the rest). 
Another possibility is that the Ca II flux in the high- and low-mass samples, while ultimately 
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related to accretion in both cases, is not probing the same physical region in the two mass 
regimes. Among the CTTs sample we have used, most of the sources (except the weakest 
accretors) evince broad Ca II emission likely arising in the magnetospheric infall (MHC98). 
Among the low-mass objects, such broad emission is seen only in GM Tau and CIDA-1; in 
the rest, the emission is narrow, and probably originates primarily in the accretion shock 
region on the (sub-)stellar surface. If so, then differences in the Ca II flux-M correlation 
might be expected between the two mass regimes. 

Further exploration of this issue is beyond the scope of this paper, but important to 
pursue in the future. For now, without a better idea of the underlying reasons for the shift 
between the low- and high-mass objects, we prefer the better fit given by equation [3] for the 
low-mass sample. However, we will conservatively assume a factor of ~5 uncertainty in the 
M inferred from this fit: this amply covers both the small observed scatter in M around the 
fit, as well as the systematic differences in M implied by equations [1] and [3] (as we will see 
in §5.3 (Table 3), the M derived from the two equations diverge by at most a factor of ~3 
for the Ca II fluxes observed in our low-mass accretors; this arises mainly from the ~0.4 dex 
difference in the intercept of the two equations, and slightly from the very small difference 
in their slopes). Similarly, for the higher- mass CTTs, we prefer the fit given by equation 
[2] over that in equation [1]. However, we do not actually use Ca II fluxes to calculate M 
for CTTs in this paper; all the CTTs accretion rates in our subsequent analysis are adopted 
from the literature. It is worth noting that these CTTs rates, derived from either veiling 
measurements or Ha modeling, also have uncertainties of a factor of 3-5 (e.g., Gullbring et 
al. 1998; MHCBH03), comparable to the errors we adopt for our low-mass sample. 

The strong correlation between Ca II flux and M in low-mass accretors enables us to 
calculate their accretion rates directly from their observed Ca II emission. We accomplish 
this in the next section. It is important to remember, however, that we do not see Ca II 
in all our accreting VLMS and BDs, only in the majority. Thus, while this technique for 
estimating M allows us to avoid laborious Ha-modeling in many cases, such modeling is 
still required in very weak accretors where Ca II emission is absent. In particular, we are 
only sensitive to line equivalent widths > O.lA. Converting this to a a flux limit (different 
for various spectral types), and using equation [3], yields a cutoff in measurable M of ~ 
1 q-io.7_ 1 q-ii.7 Mq ^-1^ going from a S p ec t ra l type of M5 to M9.5. For still lower M , Hct 

modeling is necessary (though its usefulness runs out as well, below about 1CT 12 M yr _1 ). 
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5.3. Accretion and (Sub-) Stellar Mass 

Using the above relationships, we can now calculate accretion rates for all objects in 
our low-mass sample with measured Ca II A8662 emission. The calculation is rendered 
even simpler than in CTTs, because the small accretion rates in VLMS and BDs leads to 
minimal veiling (r Ca ii), as mentioned earlier, especially in the red part of the spectrum where 
Ca II occurs (WB03). Thus we can usually assume Tcaii = Fcont x EWcaii, dropping the 
veiling term that appears in the analogous equation in §3. Of course, the absence of veiling 
around Ca II should be checked for: we have done this for our accretors with Ca II emission, 
by comparing their high-resolution spectra to those of non-accreting targets with similar 
spectral types, and indeed find negligible veiling in most cases (r Ca ii J$0.1; in IC348 382, 
MHCBH03 find rc a ii — 0.2 whereas we see none, possibly because MHCBH03 assume an 
intrinsic spectral type of M6.5, while we use the most recently determined type of M5.5). 
Veiling is present, though small, in only one of our sources - LS-RCrA-1 (~0.15) - as found 
(in the same spectrum we use) by Barrado y Navascues, Mohanty & Jayawardhana 2004. 

The accretion rates we infer from our observed Ca II fluxes in the low-mass objects, 
using both equation [1] (simultaneous fit to the CTTs and low-mass sample) and equation 
[3] (separate fit to the low-mass sources), are cited in Table 3 8 . We derive M for both the 
"bona-fide" as well as "probable" and "possible" accretors; for accretors without discernible 
Ca II emission, we quote M upper limits based on a Ca II equivalent-width detection-limit 
of O.lA. As advertised in the last section, the M from equations [1] and [3] vary by at most 
a factor of ~3 over our low-mass sample; our factor of ~5 adopted uncertainty in M covers 
both sets of values (though we prefer those from equation [3]). Independent of which fit is 
used to derive M , it is clear that accretion rates in our low-mass targets are generally far 
lower than in higher-mass CTTs. To quantify this fall-off in M with mass over the largest 
possible range in both parameters, we pool together 5 samples: (%) the objects used to derive 
the Tcaii-M relationship of equations [l]-[3], listed in Table 2; (%%) the accretors in our 
sample with M deduced from equation [3], listed in Table 3; (Hi) a few low-mass accretors 
from MHCBH03 not observed by us; (iv) additional CTTs from Gullbring et al. 1998, 



8 Note that for 2 objects - GM Tau and CIDA-1 - our Ca II fluxes have been combined with the M 
derived by WB03 to derive the M -Ca II relationships in the last section. We cannot therefore use these 
relationships to calculate new M for them from the same Ca II fluxes; for these two sources, we adopt the 
M derived by WB03 (cited in both Tables 2 and 3). For one further object - IC34 336 - we do not see any 
Ca II emission, or any veiling; however, its TiO bands do not agree very well with those of other sources with 
similar spectral types. This points to a possible small spectral typing uncertainty (of ~ 1 spectral subtype; 
our type is from Luhman 2003a); under the circumstances, we do not derive an M upper limit from the 
Ca II non-detection, but adopt instead the M found MHCBH03 from Ha modeling (cited in Table 3). 
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with M from ?7-band veiling; and (v) additional CTTs from White & Ghez 2001 (hereafter 
WG01), with M from UV excess. Objects in samples (iii)-(v) are all listed in Table 4. 

Some details about this combined sample are worth pointing out before proceeding 
further. (1) There are a few low-mass objects (IC 348 165, 205, 382 and 415), with M 
and Ca II measurements by MHCBH03, that we have independently observed as well. For 
these, we have used the MHCBH03 values in our derivation of eqn.[l,3], since MHCBH03 
derive M independently of Ca II fluxes. However, we have then applied eqn. [3] to our Ca II 
measurements to derive new M ; for these objects, it is these new M (cited in Table 3, not 
very different from those of MHCBH03 in Table 2) that we use in our M -mass analysis. 
Conversely, for MHO-5 and IC348 336, which do not show any Ca II emission in our spectra 
but do have M derived from Ha profiles by MHCBH03, we use their values. (2) All masses 
are from the Lyon tracks. WG01 already cite Lyon M.*, so these are adopted unchanged; 
MHC98, Gullbring et al. 1998, and MHCBH03 cite M* from DM tracks, so these are 
modified to Lyon values (technique described in §3). (3) All adopted M are consistent with 
Lyon tracks. M from veiling are proportional to Ai*/7Z*; thus M values from MHC98 and 
Gullbring et al. 1998, based on DM tracks, are modified to Lyon masses and radii, while M 
from WG01, based on Lyon tracks, are adopted unchanged. M derived through Ha modeling 
are rather insensitive to At* and so M from MHCBH03 are also adopted unchanged. (4) 
The WG01 sample comprises a number of binaries in Taurus. Since they resolve the binaries 
with HST, and calculate M and M* separately for each binary component, we are justified 
in including their measurements here (unlike in the MHC98 and WB03 cases mentioned in 
the last section, where we exclude the binaries because the components were not resolved, 
creating possible errors in the derived M ). However, for a number of sources, they infer 
spectral types, and thereby masses, only from photometry. Such types and masses are prone 
to error; we have excluded these sources from our analysis. (5) We have not included here 
the study by N04, which cites M inferred from near-IR spectroscopy. The reason is that, 
for their p Oph sources, spectral types and other (sub-)stellar parameters have been derived 
through IR spectral modeling; we cannot compare these quantitatively to our objects, which 
have values determined through completely different (optical) methodologies. However, most 
of the sources examined by N04 (all their Chal targets, and two of their p Oph ones) have 
been independently observed by us in the optical, and are thus already present in our sample; 
excluding N04's IR study does not therefore significantly impoverish our sample size. 

The final M -Ai* relationship is plotted in Fig. 6. For all our low-mass objects, we plot 
the M (or upper limits) computed via equation [3]. It is immediately clear that the accretion 
rate decreases sharply with mass. A linear fit to the plot, excluding all points with only upper 
limits in accretion rate, yields M oc .M* 2 , all the way from higher-mass stars to the lowest 
mass BDs. Using equation [1], instead of [3], to compute M in the low-mass sample produces 
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a nearly identical correlation (not plotted), with M oc TVf* 1 ' 85 . The small difference in the 
two relationships is not significant: as Fig. 6 shows, there is considerable scatter in the plot, 
both in the high- and low-mass regimes, which makes any attempt to precisely pin down the 
M relationship physically unjustifiable (at least with the information at hand). Under 
the circumstances, M oc M. 2 appears the simpler approximation to adopt at present. This 
is not to say that the dispersion in the data is necessarily spurious. Stochastic errors in our 
masses are unlikely to explain the scatter: while we might expect systematic uncertainties 
in our masses (since they are based on theoretical evolutionary tracks), the random errors 
should be a factor of <2 (since, regardless of any systematic offset in our masses, a well- 
determined spectral type should still correspond, from general principles, to a narrow range 
in T e ff and hence mass, over the relatively small age-span of our sample) . Similarly, in both 
our analysis and the other studies we have used, M errors are estimated to be a factor of 
~5, while the scatter in Fig. 6 is much larger: M varies by ±1.5 orders of magnitude at 
fixed mass. A significant part of the observed scatter is thus probably real, reflecting real 
variations in accretion rate for a given mass. Nevertheless, in spite of this dispersion, Fig. 6 
presents clear evidence of an overall trend of declining accretion rate with decreasing mass. 

A similar M -M.* correlation has recently been noted by other investigators as well. 
WB03 find their data to be consistent with a slope of M oc .M*; however, their study contains 
only a single object below the sub-stellar boundary with measured M . Their relationship 
holds only when BDs with M upper limits in their study are used in the fit, and thus is 
not inconsistent with our steeper slope. MHCBH03, who do have some BDs with measured 
accretion rates, find M oc M* 2 , in agreement with our results. Similarly, though N04 do not 
cite an explicit algebraic relationship, their stated numbers (M decreasing from ~10~ 8 to 
3xl0 -10 M yr _1 , going from masses around O.3-1M to <O.1M ) are also consistent with 
M oc M. 2 . However, even the latter two studies do not extend very far into the sub-stellar 
regime. In particular, the N04 investigation, which has been the most complete one prior to 
the present work, and includes their own results from infrared spectroscopy as well as data 
from WB03, MHCBH03 and various other studies, still contains only ~10 BDs (>M6.5 on 
the Luhman scale) with accretion signatures, all but one of which are M7 or earlier. Our 
present sample, which includes 7 new accreting BDs (and 2 more possibly accreting ones), 
all but one of which range from M7 to M9, not only confirms the N04 and MHCBH03's 
mass-accretion rate relationship, but shows that the correlation extends unchanged over the 
entire sub-stellar domain, down to nearly the deuterium-burning limit. 

The physical basis for this correlation, however, remains unclear. One obvious possibility 
is that the disk ionization decreases as one moves to cooler, lower-mass objects that produce 
less ionizing flux. If the magneto-rotational instability is the underlying source of disk 
viscosity driving mass accretion, as seems to be the gathering consensus, then a declining 
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disk ionization may be expected to produce lower M . This qualitative argument suggests 
that a fall-off in accretion rate with mass is perhaps not too surprising. However, it does not 
explain why the fall-off should assume the particular functional form M oc M.* . 

Alternatively, Padoan et al. (2004) have recently suggested that the accretion rate 
is fundamentally controlled by Bondi-Hoyle accretion, wherein a forming star+disk system 
gathers mass as it moves through the surrounding medium. The Bondi-Hoyle infall rate is: 

where p^, and are the gas density, sound speed and velocity of the surrounding 
material relative to the star-disk system, at large distances from the system. Ai is formally 
the mass of the star+disk; for a CTT-like object, this is almost entirely the mass of the central 
star (or BD), so that Ai ~ Ai*. If Mbh were equal to the disk accretion rate M , then 
equation [4] would immediately imply M oc At* 2 , as indicated by the observations. In reality, 
though, M B h is the rate at which surrounding material lands on the star+disk system. To 
explain the data, therefore, this scenario implicitly assumes a steady-state wherein the mass 
flux through the disk balances the mass flux falling on it at all times. Whether this can 
always be achieved remains to be seen (e.g., one might imagine instead a situation where 
excess mass builds up in the disk, and is then removed by episodic increases in the disk 
accretion rate; M would then equal Mbh only in the time-averaged sense, not at all times). 

Padoan et al. also point out that if their hypothesis is correct, then an intrinsic scatter 
can be expected as well in the M -At* relationship: this would arise from random variations 
in the stellar velocities, gas densities and sound speeds, all of which affect M by equation 
[4]. Padoan et al. then ascribe the empirical scatter seen by N04 in the M -Ai* correlation, 
and reproduced in our Fig. 6, to such variations in physical conditions. Their numerical 
simulations indicate resulting M fluctuations of ~ 1 dex, roughly consistent with the ~ 
±1.5 dex envelope seen in our sample (Fig. 6). However, this explanation of the observed 
scatter appears problematic, at least for the objects examined to date. In both our and N04's 
investigation of M versus of mass, the vast majority of the higher-mass CTTs (and many of 
the accreting VLMS and BDs as well) reside in Taurus; observationally, these objects still 
appear lined-up along the dense gas filaments in which they were presumable born (Hartmann 
2002; Briceno et al. 2002). To be sure, these sources are unlikely to still be buried within 
the filaments, since their extinctions are generally relatively low; nevertheless, they do not 
appear to have moved very far from their birthplace. Naively, then, one would expect their 
surrounding physical conditions to vary much less than over the whole cloud. However, these 
sources still exhibit ~3 orders of magnitude scatter in accretion rate at a given mass. This 
is hard to understand in Padoan et al.'s scenario: indeed, to explain an M dispersion of this 
magnitude, the latter authors invoke varying conditions over a large cloud volume certainly 
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not restricted to the vicinity of the natal filaments. Finally, denser surrounding material, 
and hence larger extinctions, might be expected towards stronger accretors in Padoan et 
al.'s scheme. This is not apparent in the Taurus sample, however. At a given mass, some of 
our CTTs differ by orders of magnitude in M while exhibiting similar Ay; similarly, most 
of the accreting Taurus BDs have Ay ~ while some non-accreting BDs in the same region 
have much higher extinctions. In summary, the relevance of Bondi-Hoyle accretion to the 
observed M -M.* relationship is not yet completely clear. 

5.4. Accretion and Age 

In Pre-Main Sequence stars, the accretion rate gradually falls off with age as the disk 
evolves and disk material is depleted (Hartmann et al. 1998). Consequently, though WTTs - 
young stars with no measurable accretion - are found even in very young star-forming regions 
(presumably reflecting a spread in disk masses and accretion rates at any given age, plus 
a spread in ages within a given region), the fraction of CTTs - stars with clear signatures 
of disk accretion - clearly declines as one moves to significantly older regions. For example, 
while the CTTs fraction in the ~1.5 Myr old Taurus region is ~50% (Kenyon & Hartmann 
1995), it is only ~15% in the ~10 Myr-old TW Hydrae Association (Muzerolle et al. 2000), 
and ~1% in the ~15 Myr-old UCL+LLC sub-groups of the Sco-Cen Association (Mamajek, 
Meyer & Liebert 2002). It is of compelling interest to ask whether a similar decline in 
accretor fraction is also apparent in the VLMS and BD regime. 

Previous studies of young low-mass objects suggest that such a decline does exist. 
JMB03 find that the proportion of VLMS and BDs with discernible accretion in high- 
resolution spectra decreases from the younger IC348 and Taurus star-forming regions to 
the older Upper Sco association. MHCBH03 find a comparable drop as well, from IC348 
and Taurus on the one hand to Upper Sco and a Ori on the other. Through a combination 
of infrared and optical spectroscopy, N04 come to a similar conclusion, finding that the ac- 
cretor fraction in the very young p Oph region significantly exceeds that in the somewhat 
older Chal cluster. Finally, applying their Ha equivalent-width criterion to low-resolution 
spectra to distinguish between accretion and activity, Barrado y Navascues & Martin (2003; 
hereafter BM03) also find a steady decline in the proportion of low-mass accretors in going 
from younger (Taurus, IC348, Chal) to older (Upper Sco, a Ori, TWHyA) regions. 

With the largest sample of VLMS and BDs observed with high-resolution spectroscopy 
to date, we can address this issue in more detail. In particular, we wish to: (1) derive the 
accretor fraction in VLMS and BDs, and identify any signs of a decline with age; and (2) 
compare the accretor fraction in this low-mass sample to that in higher-mass T Tauri stars, to 
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examine whether or not the accretor fractions in a given cluster, and the trend with age, are 
similar in the two mass ranges. The spectral type ranges used to define "low-mass objects" 
and "higher-mass T Tauri stars" for this analysis are the same as employed throughout the 
paper: >M5 (< 0.15 M Q ) for the former, and K0-M4 (2-0.25 M ) for the latter. Moreover, 
of the young clusters and star-forming regions included in this paper, we will restrict our 
attention here to p Oph, Taurus, IC348, Chal and Upper Sco; only a total of 4 objects at 
>M5 are currently known in the TW HyA and R CrA regions, too few for any meaningful 
statistics. 

Our analysis proceeds as follows. We first derive accretor fractions in the low-mass 
sample, by combining the results of our detailed accretion analysis (§5.1) with those of other 
similar recent studies. This is described in (i) below, and provides our best estimate of the 
frequency of low-mass accretors in the examined regions. However, these numbers cannot 
be compared directly to accretor fractions in higher-mass TTs, since most of the latter have 
not yet been rigorously analysed for the presence of accretion. Instead, we estimate accretor 
fractions in the TTs using the simple BM03 criteria mentioned earlier (§5.1), and compare 
these to accretor fractions in the low-mass sample re-derived using the same criteria. These 
calculations are described in (ii) and (Hi) below. While the latter estimates are less precise 
than those in (i), they do allow a valid comparison between the low and higher- mass regimes. 

(i) Low-mass (> M5) accretor fractions, from detailed accretion analysis: All our high- 
resolution spectroscopy sources (Table 1) are included in this > M5 sample, except MHO-9 
(spectral type M4.25). MHCBH03 and N04 have also conducted detailed optical and NIR 
studies of accretion; from these, we select 9 M5 and later sources in Taurus and IC348 
(MHCBH03: CIDA-14, FN Tau, MHO-7, MHO-8, V410 Anon 13, V410 X-ray 5a, IC348 
173, 407, 454) and 8 in p Oph (N04: ISO-Oph 23, 32, 33, 102, 160, 176, 193, GY 10) 
that we have not observed 9 . GY 310 and IC348 256, which we designate as "probably" 
accreting in §5.1, are counted here as accretors. However, we conservatively consider the 



9 We have not included N04 sources in our previous M vs. M* analysis, since N04 derive A4* and 
M through NIR calibrations of T e // , luminosity and extinction, which are not directly comparable to our 
optically calibrated values without analysis beyond the scope of this paper. However, in the present analysis, 
we are interested only in whether the spectral types of the N04 sources are >M5 or not, and whether they 
are accreting ot not. Through NIR comparisons, N04 find all their p Oph sources to be >M6. In the present 
study, we use optical spectral types; Luhman & Rieke 1999 have shown that optical types for young objects 
are generally ~1 sub-type earlier than NIR types for classes M6-M7, and <~l-2 earlier for classes >M8. 
Using this conversion, the N04 p Oph sources would still have (optical) spectral types >M5; this expectation 
is borne out by comparing N04's types to those of Luhman & Rieke for 4 sources in common. Moreover, 
the qualitative identification of accretion by N04 in their sources is independent of the precise spectral type. 
Thus, their sources can be included in our present analysis of accretor fraction at spectral types >M5. 
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"possible" accretors KPNO-Tau 4, KPNO-Tau 11 and CFHT-BD Tau 4, as well as the 
possible spectroscopic binaries KPNO-Tau 14 and DENIS 161929, as non-accretors in the 
present analysis. For the MHCBH03 and N04 sources, we adopt the accretor / non-accretor 
status found by those authors. Finally, since IC348 and Chal appear to be at roughly the 
same age of ~2 Myr, we do not differentiate between them in this accretion-age analysis, 
but instead combine their members into one group. This serves three purposes: it provides 
a single large sample at ~2 Myr, dilutes our bias towards accretors in the M5-5.5 range in 
IC348 (see §2), and offsets the relatively small number of BDs (type > M6) observed in Chal 
alone. Accretor fractions are calculated as N acc /N tot) and errors as \fN acc / 'N tot . 

The resulting accretor fractions are: 58±22% (7/12) in p Oph, 32±10% (10/31) in 
Taurus, 32±10% (10/31) in IC348 + Chal, and 5±5% (1/20) in Upper Sco. These are 
plotted in Fig. 7 (black triangles). We see a distinct decrease in the frequency of accretion, 
moving from p Oph at <1 Myr, to Taurus, IC348 and Chal at 1.5-2 Myr, to Upper Sco at 5 
Myr. Given the error bars, the difference between p Oph on the one hand and Taurus, IC348 
and Chal on the other is not definitive: comparable accretor fractions in these regions is not 
very surprising, since they are all close in age at <2 Myr. However, a large fall-off in accretor 
frequency by Upper Sco ages is clear: accretion appears to have generally halted by ~5 Myr 
(or more accurately, dropped below measurable levels by this age, corresponding to M < 
10~ 12 M yr -1 given our Ha 10% width > 200 km s" 1 criterion for identifying acretors). 

(ii) Low-mass (> M5) accretor fractions, from BM03 criteria: BM03 define spectral 
type-dependent Ha equivalent width criteria to identify accretion. We apply this technique 
to the same low-mass sample as above, to re-derive accretor fraction estimates that can be 
compared to T Tauri stars. However, no Ha measurements are available for 6 of N04's p Oph 
sources included in our above sample; these are dropped from the present analysis. Moreover, 
the BM03 criteria are defined for, and thus best applied to, low-resolution spectra. We have 
thus used Ha low-resolution equivalent widths collected from the literature for most of our 
low-mass sources 10 . For a handful of objects, no low- resolution Ha data are available; for 
these, we were forced to apply the BM03 scheme to high-resolution Ha equivalent widths 
(GY 5, 37, 310 and USco 130, 131 from Table 1; ISO-Oph 32 and 102 from N04). 

The BM03 method implies low-mass accretor fractions of: 67±33% (4/6) in p Oph, 
45±12% (14/31) in Taurus, 52±13% (16/31) in IC348 + Chal, and 20±10% (4/20) in Upper 



10 pOph: Luhman & Ricke 1999 (GY 141). Taurus: Bricefio et al. 2002 & Luhman ct al. 2003b (KPNO 
sources); Martin ct al. 2001 (CFHT sources); Bricefio et al. 1999 (CIDA-14); White et al. 1999 (GG Tau 
Ba,b); Bricefio et al. 1998 (MHO sources, V410 Anon 13, V410 X-ray 5a); Kenyon ct al. 1998 (CIDA-1, FN 
Tau); Strom & Strom 1994 (V410 X-ray 3). 7C348: Luhman et al. 2003a. Chal: Luhman 2004. UpperSco: 
Ardila, Martin & Basri 2000 (USco sources); Martin, Delfosse & Geieu 2004 (DENIS sources). 
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Sco. These values are plotted in Fig. 7 (red circles), and are systematically higher than those 
inferred in (i) from a detailed study of accretion. Compared to our high- resolution analysis, 
the BM03 scheme is more prone to mis-classifying high activity as accretion; conversely, 
it may be better at identifying high-inclination, low-M accretors, whose high-resolution 
Ha profiles would suggest activity alone (as discussed in §5.1). Both effects yield a less 
conservative (higher) estimate of accretor fraction from BM03 criteria than derived in (i). 
Nevertheless, the two sets of estimates are comparable within the errors, and the BM03 
fractions reveal the same trend with age: roughly constant accretor frequencies in low-mass 
objects from <1 to 2 Myr, followed by a large drop by the 5 Myr age of Upper Sco. 

(Hi) Higher- Mass T Tauri (K0-M4) accretor fractions, from BM03 criteria: Finally, 
we use the BM03 technique to derive accretor fractions for higher-mass T Tauri stars, for 
comparison to the low-mass sample above. BM03 themselves list accretor frequencies for T 
Tauri stars in various regions. However, we prefer not to use their values, for three important 
reasons. First, their 'T Tauri' sample extends down to spectral type M5.5, which overlaps 
with our low-mass sample of >M5; as such, it does not provide a clean comparison between 
objects near the substellar boundary and those that are significantly higher in mass. Second, 
their p Oph T Tauri sample includes sources not just in the dark core of this region (where our 
low-mass p Oph sources reside), but also in the surrounding areas. Martin et al. (1998) have 
shown that the ratio of CTTs to WTTs decreases sharply outside the dark core; indeed, this 
probably accounts for the puzzlingly low CTTs fraction that BM03 find, and draw attention 
to, in p Oph as a whole. For a valid comparison between the low and high-mass regimes, 
therefore, we must select spatially proximate objects in p Oph. Finally, the BM03 'Upper 
Sco' sample includes sources in the Sco-Cen UCL+LLC sub-groups as well; since the latter 
are expected to be significantly older than the Upper Sco region alone (Mamajek, Meyer & 
Liebert 2002), they are unsuitable for inclusion in our age analysis here. 

We have therefore compiled anew a list of K0-M4 stars in the pertinent star-forming 
regions and young clusters. In p Oph, all our low-mass objects lie within a ~35'x35' area 
of the dark core originally studied by Casanova et al. 1995. We select K0-M4 stars in this 
region from Bouvier & Appenzeller 1992 and Martin et al. 1998. For Taurus, we select stars 
from Briceno et al. 2002 and Luhman et al. 2003b (since most of our Taurus low-mass 
objects lie in the Taurus fields analysed in these two studies). For IC348 and Chal, we select 
stars from Luhman et al. 2003a and Luhman 2004 respectively; the latter studies provide 
the currently most comprehensive list of stars in these regions. Finally, for Upper Sco we 
select stars from the studies by Walter et al. 1994, Preibisch et al. 1998 and Prebisch et al. 
2002, which together provide the most complete list of stars with available Ha data in this 
cluster to date. We have then compiled Ha equivalent widths for all these stars from various 
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sources in the literature 11 - mostly low-resolution, but high-resolution in a handful of cases 
where low-resolution data is unavailable - and applied the BM03 accretion criteria. 

The resulting K0-M4 T Tauri accretor fractions (i.e., CTTs fractions) are: 50±16% 
(10/20) in p Oph; 59±9% (42/71) in Taurus; 33±6% (29/87) in IC348 and 44±8% (28/63) 
in Chal (so 38±5% in IC348 + Chal); and 7±2% (12/170) in Upper Sco. These are plotted 
in Fig. 7 (blue squares), and immediately point to two results. 

First, the decline in CTTs fraction is similar to that in the low-mass sample: while the 
CTTs frequencies are comparable in p Oph, Taurus, IC348 and Chal, they fall off steeply by 
the age of Upper Sco. A closer look reveals that our CTTs fraction in ChaI+IC348, while 
consistent with that in p Oph, is in fact somewhat lower than in Taurus. The significance of 
this is unclear with the present errors; better number statistics, especially in p Oph, can tell 
us whether there is a decline by the age of IC348 and Chal, compared to the younger p Oph 
and Taurus, or whether Taurus has a larger CTTs fraction relative to both the younger p Oph 
and older IC48+ChaI regions. At any rate, the data at hand show very clearly that the CTTs 
fraction has declined vastly by the age of Upper Sco, in agreement with the trend observed 
in the low-mass sample. Second, in any given region, the high- and low-mass regimes evince 
comparable accretor fractions, within the errors (using BM03 criteria for both) 12 . 

In summary, VLMS and BDs appear quite similar to higher mass T Tauri stars, both in 
the frequency of accretion at a given age, and in the overall accretion timescale (of order ~5 
Myrs) . It is worth noting that a few CTTs are also known in the TW Hydrae Association, 
which at 10 Myrs is older than Upper Sco. We have not analysed TW HyA above, since only 
three low-mass objects have been identified in it so far (all BDs; Table 1); nevertheless, at 
least one of these sources is accreting as well (2MASS 1207-3932; see §5.1). This supports 
our conclusion that accretion in low-mass objects can last just as long as in T Tauri stars. 



5.5. Accretion and Disks 

Lastly, we examine the relationship between disk-accretion signatures, and infrared 
excesses directly indicating the presence of dusty disks. Jayawardhana et al. (2003; hereafter 



u pOph: Bouvier & Appcnzcller 1992; Martin et al. 1998. Taurus: Hartmann et al. 1991; Hartigan, 
Strom & Strom 1994; Martin et al. 1994; Mathieu 1994; Kenyon et al. 1998; Monin, Menard & Duchene 
1998; Briceno ct al. 1998; Bricerio et al. 1999; White & Hillenbrand 2004. IC348: Luhman et al. 2003a. 
Chal: Luhman 2004. UpperSco: Walter et al. 1994; Preibisch et al. 1998; Preibisch et al. 2002. 

12 In Upper Sco, there appears to be a slight difference in accretor frequency between the two mass regimes, 
using the BM03 scheme, but the difference is not significant with our present low-mass number statistics. 
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JASH03) have carried out an L'-band survey of various young clusters and star-forming 
regions; they (conservatively) infer the presence of a disk when the K — V excess is > 0.2 
mag. After making a couple of minor revisions to their results (to account for updated 
spectral types and membership status for a few objects), we find that the frequency of disks 
in M5 and later objects is 67±27% (6/9) in p Oph, 55±25% (5/9) in Taurus, 53±17% (10/19) 
in IC348 + Chal, and 60±35% (3/5) in Upper Sco. These values are consistent, within the 
errors, with a roughly constant disk fraction of 50-60% from p Oph down to Upper Sco ages. 
This clearly differs from the trend in accretor fraction in the same regions, which, as we 
showed above, distinctly decreases by the 5-Myr age of Upper Sco. 

To explore this issue further, we compare individual objects in common between JASH03's 
sample and ours. First, the 6 accretors we have in common in p Oph, Taurus, IC348 and 
Chal all evince K — V excesses > 0.2 mag. Only one accretor in our study - 2MASS 1207- 
3932 in TWHyA - does not exhibit any K — V excess. However, Apai et al. (2004) have 
recently detected disk-related mid-infrared emission from this source. Their data suggest 
significant grain-evolution in the disk, which might account for the lack of an L'-band ex- 
cess; alternatively, L'-band emission may be suppressed by a inner hole in the disk. In any 
event, there is no evidence for accretion signatures arising in the complete absence of any 
disk excess, a situation that would be quite difficult to explain. 

However, there are a number of sources with K — L' excess that do not have any 
measurable accretion (see Table 1). For example, the Chal sources ChaHa 1, 7, 8, 9 and 12, 
and the Upper Sco objects USco 104, 112 and 128, all have excesses greater than 0.2 mag, and 
are thus likely to possess disks; none of them appears to be accreting in our high-resolution 
spectra. This strongly implies that significant amounts of disk material can persist around 
young VLMS and BDs, even after accretion has halted, or at least declined to unmeasurable 
levels. This is consistent with our statistical argument above, showing that the disk fraction 
remains high even when the accretor fraction has diminished substantially with age. This 
situation in the low-mass sample is analogous to the transition from the CTT to WTT phase 
in higher-mass stars: WTTs by definition do not have measurable accretion, yet often evince 
disk excesses (e.g., Haisch, Lada & Lada 2001). 

6. Concluding Remarks 

Our results may be summarized thus: (1) Classical T Tauri-like accretion is widespread 
in the very low-mass stellar and sub-stellar domain, extending to the least massive brown 
dwarfs examined so far; (2) accretion in these objects can be accompanied by mass outflows, 
as in CTTs; (3) the accretion rates decline sharply with mass, approximately as M oc Ai* 2 
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(with considerable scatter); this functional form for the decline extends not only to the lowest 
mass stars, but also over the entire sub-stellar regime; (4) the accretor frequency decreases 
with age in very low-mass stars and brown dwarfs, as in higher-mass CTTs; (5) however, 
like in CTTs, accretion can continue in the low-mass regime for up to 10 Myrs; at any given 
age the accretor fraction in both mass regimes is comparable, implying similar accretion 
timescales; and (6) disks around objects near and below the hydrogen-burning limit persist 
beyond the main accretion phase, mirroring the transition from the classical to post-T Tauri 
phase in solar-type stars. These phenomenological similarities suggest that very low-mass 
stars and brown dwarfs form in the same way as higher-mass stars do. We find no evidence to 
support a change in formation mechanism between the two mass regimes, such as postulated 
in the 'ejection' scenario. 

Nevertheless, it is worth remembering that the optical accretion diagnostics explored 
here, as well as the infrared excess measurements that attest to the ubiquity of sub-stellar 
disks, probe only the innermost disk radii. It is still possible that outer disks are indeed 
truncated through dynamical ejections, while the inner regions survive for a long time merely 
by virtue of the small accretion rates in the low-mass regime. Estimates of disk sizes (and 
masses) in brown dwarfs are required for a more conclusive resolution of this question. The 
detection of optically thin sub-millimeter and millimeter emission, which would indicate the 
existence of cool outer disks, and also provide a direct estimate of dust (and, by extrapolation, 
disk) masses, would be vastly helpful in this regard. Such emission has already been detected 
in two young brown dwarfs (Klein et al. 2003); larger surveys are currently underway, and 
should greatly clarify this issue of origins in the near future. 
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Fig. 1. — Ha profiles for all our accretors, including the 5 objects designated "probable" 
and "possible" accretors in the text (§5.1). Dashed horizontal lines indicate the continuum 
and the Ha 10% width. The emission is usually strong with broad line-wings, and often 
asymmetric. Note that eight objects with very high Ha equivalent widths are plotted on a 
y-axis scale larger than in the others; these panels are marked by a thick right-hand y-axis 
(in seven of these, the y-scale is twice that in the others; in KPNO-Tau 12, with the highest 
equivalent width in our sample, the scale is four times larger). For LS-RCrA-1 (third panel 
from end), we have obtained 3 consecutive spectra, and show here the averaged Ha profile 
(which does not vary much between the observations). For 2MASS 1207-3932 (last panel), 
we also have 3 consecutive spectra; since its Ha profile varies significantly over these ~1 hour 
time-scales, we overplot all three (thick solid, thin solid, and dashed lines). The separated 
profiles, and their discussion, can be found in Mohanty, Jayawardhana & Navascues (2003). 
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Fig. 2. — Ha profiles for chromospherically active, non-accreting sources in our sample (only 
new objects are shown; the rest have been published previously - see §5.1 and Table 1). 
Compared to the accretors (Fig.l), the emission here is usually much weaker, narrower, and 
symmetric. Four objects with extremely small emission are plotted on a y-axis scale half that 
in most others; these panels are marked with a dashed right-hand y-axis. Two sources with 
very strong (but still narrow and symmetric) emission are plotted on a y-scale larger than the 
others: by a factor of five in Cha Ha 1 and a factor of two in Cha Ha 7 (panels marked with 
a thick right-hand y-axis). The Ha profile in KPNO-Tau 8 is comparatively broad, due to its 
high v sim (45 km s" 1 ); KPNO-Tau 14 and USco DENIS 161929 appear to be spectroscopic 
binaries from our cross-correlation analysis, probably accounting for the somewhat broad 
profile in the former and asymmetry in the latter. Emission is also asymmetric in Cha Ha 
11 and 12; they might be accreting, but do not evince any other signatures of doing so in 
our data (§5.1), so we classify them as active for now. Note that the emission at ~ 6569A m 
some panels is a sky-line, not stellar, that remained unsubtracted. 
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Fig. 3. — Observed Ca II A8662.14 emission in our accretors. The emitting objects are 
plotted in black, the emission line is marked with an arrow; the horizontal dashed line is 
the normalized continuum. For comparison, we overplot in grey the same spectral region for 
KPNO-Tau-5, a young, non-accreting M7.5 object in Taurus with good S/N and no Ca II 
emission. The comparison shows that most of the bumps and wiggles in our accretors, in 
the continuum around the Ca II line, are real (there are some mismatches due to varying 
spectral type - our accretors range from M5 to M9.5 - and low S/N in a few sources). 
In particular, the excellent agreement between KPNO-Tau 5 and our weakest emitters - 
KPNO-Tau 6, KPNO-Tau 7 and CFHT-BD-Tau 4 - at all points in the continuum except 
only at the Ca II position, shows that Ca II emission is indeed present in the latter objects. 
Note that the dip in the spectrum of KPNO-Tau-5 at the approximate position of Ca II is 
not due to photospheric absorption in this line: the photosphere is much too cool for this. 
Instead, this dip is just another wiggle in the continuum that happens to overlap with the 
Ca II position. The observed emission in most of our accretors is fairly narrow, with FWHM 
< 50km s _1 . Nevertheless, as described in §5.1, the emission is seen only in the accretors; 
in merely chromospherically active low-mass sources, the spectrum simply looks like that of 
KPNO-Tau 5, which is active as well (evincing chromospheric Ha). 
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Fig. 4.— Forbidden [O I ] A6300 emission in the accreting M8 TW HyA brown dwarf 2MASS 
1207-3932 (thick line). The horizontal dashed line is the continuum (normalized over a much 
larger wavelength range than plotted). For comparison, we overplot the spectrum of 2MASS 
1139-3159 (thin line), which is a non-accretor with the same spectral type, also in TW HyA. 
Excess emission is clearly visible in the accretor, indicating that mass outflows can persist, 
like accretion, in low-mass BDs for up to 10 Myr. 
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Fig. 5. — Observed dependence of Ca II emission-flux on M (where M is calculated inde- 
pendently of the Ca II flux, from either veiling or Ha-modeling; §5.2). Higher-mass CTTs 
are shown as filled circles, and low-mass accretors (VLMS and BDs) as filled squares. The 
Ca II flux declines sharply with decreasing accretion rate over the entire sample, with the 
correlation extending over ~4 orders of magnitude in both flux and M . Top panel: Single 
fit to both the CTTs and low-mass accretors (equation [1], §5.2). Notice that the low-mass 
sample lies slightly above the overall fit, while most of the CTTs lie slightly below. Bottom 
panel: Separate fits to the CTTs and low-mass accretors (equations [2] and [3], §5.2). The 
sytematic offset between the CTTs and low-mass objects is now obvious, with the CTTs 
evincing lower Ca II flux at a given accretion rate. The scatter around the individual fits is 
much reduced, compared to the single fit in the top panel. See text (§5.2). 
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Fig. 6. — Derived relationship between M and mass. Low- mass accretors with M derived 
in this paper, from Ca II fluxes using eqn.[3], are shown as either filled circles (bona-fide 
accretors with detected Ca II ), bulls-eyes ("probable" or "possible" accretors with detectd 
Ca II ) or downward arrows (accretors without detected Ca II emission: M upper limits 
derived from Ca II detection limit of O.lA). Crosses are CTTs from Gullbring et al. 1998; 
pluses are CTTs in binaries from WG01; squares are low-mass sources from WB03; asterisks 
are low-mass sources from MHCBH03 (see §5.3). The thick solid black line is our formal fit 
to the data (excluding the upper limit points): M oc M.* . The two dashed lines are this 
fit vertically offset by ±1.5 dex, to denote the upper and lower envelopes of the trend. All 
sources lie within (or very close to) this envelope, suggesting that the data is consistent with 
scatter around a single functional relationship between M and Ai* such as we have derived 
(however, the data is not sufficient to prove this). The main implication is that M decreases 
sharply with mass, following the same slope from CTTs to VLMS to the lowest mass BDs: 
over nearly 2.5 orders of magnitude in mass and 4 orders in M . See text (§5.3). 
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Fig. 7. — Accretor fraction (%) as a function of age. The regions selected to represent various 
ages are marked: p Oph (<1 Myr), Taurus (1.5 Myr), IC348+ChaI (2 Myr; members of these 
two coeval regions are grouped together), and Upper Sco (5 Myr). Black triangles: accretor 
fractions in low- mass sources (spectral type >M5), from a detailed analysis of accretion 
signatures. Red circles: accretor fractions in the low-mass sources, now using BM03 criteria 
instead. Blue squares: accretor fractions in higher mass stars (K0-M4), again using BM03 
criteria. At a given age, the three symbols are slightly offset horizontally for clarity. We see 
that (1) accretor fractions from BM03 criteria are systematically higher than from detailed 
analysis, in the low-mass sample; (2) the accretor fractions in both the low-mass sample 
(using either BM03 criteria or detailed analysis) and the higher-mass stars fall sharply by 5 
Myr; and (3) at any fixed age, the accretor fractions in the low- and high-mass regimes are 
comparable (using BM03 criteria for both). See text (§5.4). 
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Table 1 — Continued 
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10.0 


15. 


208. 




1: 


0.2 


Y 




IC348 


luo48 405 


M8.0 


45: 


1.7: 


141: 








N 




IC348 


10o48 415 


M6.5 


0.0 


80. 


213. 




3.4 


1.7 


Y 




IC348 


luo48 478 


M6.25 


20.0 


3.1 


150. 








N 




IC348 


Una Ha 1 


JYu . / 5 


8.0 


173. 


151. 








IN 


disk 


Una 1 


Una Ha 2 


i\ /r p- oc 

IVlO.ZD 


20.0 


63.3 


317. 








Y 


disk 


una 1 


Cha Ha 3 


M5.5 


22.0 


24.0 


136. 








N 




Cha I 


pL„ TT _ /I 

Una Ha 4 


M5.5 


20.0 


8.3 


127. 








IN 




una 1 


Una Ha o 


M5.5 


16.0 


17.8 


111. 








TVT 

IN 




una 1 


Una Ha o 


MO. ( O 


c: a 


oy.b 


OQO 

zoz. 








v 
i 


disk 


una 1 


Cha Ha 7 


M7.75 


11.0 


52.2 


114. 








N 


disk 


Cha I 


Cha Ha 8 


M5.75 


10.0 


8.4 


106. 








N 




Cha I 


Cha Ha 9 


M5.5 


9.0 


2.4 


91. 








N 


disk 


Cha I 


Cha Ha 10 


M6.25 


10.0 


2.5 


83. 








N 




Cha I 


Cha Ha 11 


M7.25 


16.0 


12.2 


136. 








N 




Cha I 


Cha Ha 12 


M6.5 


22.5 


6.8 


107. 








N 


disk 


Cha I 


Cha Ha 13 


M5.5 


14.0 


14.3 


122. 








N 




Cha I 


LS-RCrA-1 


M6.5 


18.0 


53. 


316. 




Y 


9.6 


Y 




R CrA 


USco 40 


M5.0 


37.5 


7.9 


137. 








N 




Upper Sco 


USco 53 


M5.0 


45.0 


17.8 


175. 








N 




Upper Sco 


USco 55 


M5.5 


12.0 


7.3 


114. 








N 




Upper Sco 


USco 66 


M6.0 


27.5 


6.5 


115. 








N 




Upper Sco 


USco 67 


M5.5 


18.0 


12.9 


139. 








N 




Upper Sco 


USco 75 


M6.0 


63.0 


8.9 


212. 








N 




Upper Sco 


USco 100 


M7.0 


50.0 


13.1 


184. 








N 




Upper Sco 



Table 1 — Continued 





SpT 


i) sin % 


Ha 


Wrv 1 0% width 

1JU J-O/0 W1LIU11 


jnc± auu i o 






cx\jKjL L LUl ; 


/ • / i r 1 1 i V i f\y\\- O 
l_,U 11111 


region 






[km s _1 ] 


rAi 


[km s _1 ] 


rAi 


fAl 


fAl 


rv/Ni 






USco 104 


M5.0 


16.0 


9.4 


109. 


_ 


_ 


_ 


N 


disk 


Upper Sco 


USco 109 


M6.0 


6.0 


6.3 


84. 


- 


- 


- 


N 


- 


Upper Sco 


USco 112 


M5.5 


8.0 


9.5 


111. 








N 


disk 


Upper Sco 


USco 128 


M7.0 


0.0 


15.9 


121. 


0.3 






N 


disk 


Upper Sco 


USco 130 


M7.0 


14.0 


8.4 


111. 








N 




Upper Sco 


USco 131 


M7.0 


16.0 


14.2 


105. 








N 




Upper Sco 


USco DENIS 160603 


M7.5 


9.0 


70. 


306. 








Y 




Upper Sco 


USco DENIS 161006 


M8.5 


7.0 


13.9 


110. 








N 




Upper Sco 


USco DENIS 161103 


M9.0 


0.0 


16.9 


90. 








N 




Upper Sco 


USco DENIS 161452 


M9.0 


0.0 


5.3 


79. 








N 




Upper Sco 


USco DENIS 161916 


M8.0 


5.0 


10.3 


93. 








N 




Upper Sco 


USco DENIS 161929 


M8.0 




13.6 


135. 








N 


SB2 


Upper Sco 


USco DENIS 162041 


M7.5 


16.0 


7.4 


81. 








N 




Upper Sco 


2MASS 1207-3932 


M8.0 


13.0 


27.7 


204. 


0.8 






Y 


disk 


TW HyA 


2MASS 1139-3159 


M8.0 


25.0 


7.3 


111. 








N 




TW HyA 


TWA 5B 


M8.5 


16.0 


5.1 


162. 








N 




TW HyA 



a Accreting or not, as discussed §5.1. Non-accretors marked with 'N', "bon-fide" accretors marked with 'Y' (including MHO-5; sec [d] below), and "probable" 
and "possible" accretors marked as such (see §5.1). 

^Objects with known direct evidence for disks, from excesses in the mid-IR (2MASS 1207-3932) and/or in the near-IR L'-band (all others). Note that 
sources with no information may still have disks. 

c Our Keck spectra for CFHT-BD-Tau 2 and 4 do not cover the He I A6678 line. However, we do not see any overt evidence for He I emission in their 
published low-resolution spectra (Martin et al. 2001). 

d MHO-5 looks like a non-accretor in our spectrum, but shows forbidden [O I ] A6300, strongly indicative of an outflow, and hence accretion, in MHCMH03's 
data ([O I ] is not covered in our spectrum); we thus assume it is a bona-fide accretor (see §5.1). 
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Table 2. Objects Used to Derive Ca II Flux - Accretion Rate Relationship 



name 


op 1 


L eff 


a 

rcall 




■TCall 


T? c 
/V* 




lVl 


acc-method 


rei 






[K] 




[A] 


[log(crg s~ 1 cm ) 


1 [Rq] 


[Mq] 


[log{M Q yr- 1 )] 






AA Tau 


K7 


4050 


0.0 


0.5 


5.74 


2.14 


0.78 


-8.56 


ex 


i 


BP Tau 


K7 


4050 


0.6 


7.8 


7.13 


2.14 


0.78 


-7.71 


ex 


i 


CW Tau 


K3 


4750 


0.8 


11.3 


7.68 


2.94 


1.60 


-7.21 


ex 


i 


DG Tau 


K7/M0 


3950 


3.6 


30.3 


8.16 


2.13 


0.75 


-6.93 


ex 


i 


DN Tau 


MO 


3850 


0.0 


0.5 


5.67 


2.10 


0.70 


-8.72 


ex 


i 


DR Tau 


K7 


4050 


5.1 


21.8 


8.16 


2.14 


0.78 


-6.66 


ex 


i 


DS Tau 


K5 


4350 


0.3 


3.3 


6.84 


2.32 


1.00 


-7.72 


ex 


i 


GK Tau 


K7 


4050 


0.0 


0.5 


5.74 


2.14 


0.78 


-8.42 


ex 


i 


CIDA-1 


M5.5 


3050 


0.69 


23.6 


7.12 




0.12 


-8.5 


ex 


2 


GM Tau 


M6.5 


2950 


1.04 


14.5 


6.93 




0.073 


-8.6 


ex 


2 


MHO-5 


M6.0 


3000 


0.0 


0.1 


4.48 




0.10 


-10.8 


Ha 


3 


IC348 165 


M5.25 


3100 


0.0 


0.4 


5.16 




0.14 


-10.0 


Ha 


3 


IC348 173 


M5.75 


3025 


0.0 


0.7 


5.35 




0.11 


-10.0 


Ha 


3 


IC348 205 


M6.0 


3000 


0.1 


0.9 


5.47 




0.10 


-10.0 


Ha 


3 


IC348 382 


M5.5 


3050 


0.2 


0.3 


5.07 




0.11 


-10.8 


Ha 


3 


IC348 415 


M6.5 


2950 


0.1 


2.6 


5.90 




0.07 


-9.3 


Ha 


3 



a Veiling in the vicinity of Ca II A8662 (adopting the values calculated at 8600A for sources in rcf. (1), 8400A for sources in ref. (2), and 
8900 A for sources in ref. (3)). 

b Ca II A8662 emission flux, calculated by us by combining the observed equivalent width and veiling with the continuum flux estimated 
for the given spectral type (§2). 

c Mass and radius from the Lyon tracks, for the given spectral types (see §2). Radii are calculated only for those sources where they are 
needed to modify the M (sec [d], [c], [f] below). 

d Accretion rate from various sources, modified by us from the original values in some cases (§5.1; see notes [e] and [f] below). 

e Method used to derive M . 'ex': from optical veiling excess; 'Ha': from Ha line-profile modeling. 

f References for the spectral types, veiling, Ca II A8662 equivalent widths and M . 1: MHC98; 2: WB03; 3: MHCBH03. MHC98 use DM 
tracks to infer M from veiling, so we have modified their original M to Lyon masses and radii (§5.1); these modified M (and new masses 
and radii) are cited here. The MCHBH03 M from Ha-modcling arc adopted unchanged; we have only changed the masses they cite to Lyon 
values. The WB03 values for M (from UV excess) and mass (from Lyon tracks) are both adopted unchanged. Finally, WB03 do not supply 
Ca II A8662 equivalent widths for GM Tau and CIDA-1, but we have calculated these from the same spectra they used. 
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Table 3. Derived Masses and Accretion Rates in Our Sample 



name 


op 1 


L eff 


rcall 


■TCall 


A A 


1V1 


1V1 






[ K J 




[log(erg s^ 1 cm -2 )] 


l M ©J 


l'og(M Q yr )J 


[tog(WlQ yr JJ 


GY 5 


M5.5 


3050 


0.0 


<4.52 


0.12 


<-10.6 


<-10.8 (-10.8, -10.1) c 


GY 310 


M7.0 


2900 


0.0 


4.70 


0.06 


-10.4: 


-10.6: ( — , -9.3) c 


KPNO-Tau 3 


M6.0 


3000 


0.0 


5.18 


0.08 


-9.9 


-10.2 


KPNO-Tau 4 


M9.5 


2300 


0.0 


4.08 


0.015 


-11.1: 


-11.2: 


KPNO-Tau 6 


M8.5 


2550 


0.0 


4.34 


0.025 


-10.8 


-11.0 


KPNO-Tau 7 


M8.25 


2650 


0.0 


4.18 


0.03 


-11.0 


-11.2 


KPNO-Tau 10 


M5.0 


3125 


0.0 


<4.58 


0.16 


<-10.5 


<-10.7 


KPNO-Tau 11 


M5.5 


3050 


0.0 


<4.52 


0.13 


<-10.6: 


<-10.8: 


KPNO-Tau 12 


M9.0 


2400 


0.0 


4.62 


0.02 


-10.5 


-10.7 


CFHT-BD-Tau 4 


M7.0 


2900 


0.0 


4.40 


0.06 


-10.7: 


-10.9: 


MHO-5 


M6.0 


3000 


0.0 


<4.48 


0.10 


<-10.7 


<-10.9 (-10.8) c 


IC348 165 


M5.25 


3100 


0.0 


5.60 


0.14 


-9.5 


-9.9 


IC348 205 


M6.0 


3000 


0.0 


5.08 


0.10 


-10.0 


-10.3 


IC348 256 


M5.75 


3025 


0.0 


<4.51 


0.12 


<-10.6: 


<-10.8: 


IC348 291 


M7.25 


2850 


0.0 


<4.36 


0.05 


<-10.8 


<-11.0 


IC348 382 


M5.5 


3050 


0.0 


4.82 


0.11 


-10.3 


-10.6 


IC348 415 


M6.5 


2950 


0.0 


5.68 


0.07 


-9.4 


-9.8 


Cha Ha 2 


M5.25 


3100 


0.0 


<4.56 


0.14 


<-10.6 


<-10.8 (-10.0, -10.8) c 


Cha Ha 6 


M5.75 


3025 


0.0 


<4.51 


0.12 


<-10.6 


<-10.8 (-10.5, <-10.8) c 


LS-RCrA-1 


M6.5 


2950 


0.15 


6.50 


0.08 


-8.5 


-9.0 


USco DENIS 160603 


M7.5 


2800 


0.0 


<4.32 


0.04 


<-10.8 


<-11.0 


2MASS 1207-3932 


M8.0 


2700 


0.0 


<4.20 


0.035 


<-10.9 


<-ll.l 


CIDA-1 


M5.5 


3050 


0.69 


7.12 


0.12 




- ("8.5) d 


GM Tau 


M6.5 


2950 


1.04 


6.93 


0.073 




- (-8.6) d 


IC348 336 


M5.5 


3050 






0.10 




- (-10.0) d 



a Accrotion rates for our low-mass sample, inferred from the Ca II A8662 fluxes via equation [1] (i.e., using a simultaneous fit 
to the higher-mass CTTs and low-mass accretors; §5.2). For sources with no detected Ca II emission, upper limits for M are 
given, based on our detection limit of 0.1 A for the Ca II equivalent width. For the 5 sources classified as "probable" or "possible" 
accretors (Table 1 and §5.1), our derived M are marked with a colon. 

b Same as previous column, but now M calculated using equation [3] (i.e., using the fit to only low-mass accretors; §5.2). 

c For objects in which we have derived M upper limits, but M measurements have been made by other investigators as well, we 
cite the latter values in parantheses in addition to our upper limits. For MHO-5, we cite M from MHCBH03 (from Ha modeling); 
for GY 5, GY 310, Cha Ha 5 and Cha Ha 6, we cite M from N04 (from both Ha modeling and NIR Pa/3 luminosities, in that 
order). While MHCMH03's and N04's values are higher than our upper limits in some cases, we do not consider this a serious 
problem: the discrepancies are generally quite small, and within an order of magnitude in all but one case; they are within the 
range expected from the differences in our techniques for determining the (sub-)stellar parameters and M , as well as from possible 
M variations in time (indeed, similar differences are present even within N04's values, derived using two different techniques). 

d For GM Tau, CIDA-1 and IC348 336, we do not calculate M from Ca II fluxes, but adopt instead the M values from WB03 
(GM Tau, CIDA-1) and MHCBH03 (IC348 336), cited in parantheses. For GM Tau and CIDA-1, we have used our Ca II fluxes, 
combined with the WB03 M , to derive the M — Ca II relationships (see Table 2), so we cannot use these relationships to derive M 
again from the same Ca II fluxes. In IC348 336, we do not see any Ca II A8662 emission or any clear signs of veiling which might 
hide the line. However, its TiO bands around 8440A do not match too well those of other sources at the same spectral type. The 
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mismatch docs not seem to result from veiling, though MHCBH03 find a small possible veiling of ~0.1±0.1 at 7100 and 8900A. 
The TiO discrepancy may be due to uncertainty in its spectral type (within 1 subclass) which we have adopted from Luhman et 
al. 2003a. As such, we choose not to derive an M upper limit based on our Ca II non-detection, but adopt instead the M found 
by MHCBH03 from Ha modeling. 
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Table 4. Masses and Accretion Rates in Other CTTs and Low-Mass Accretors 



name 


SpT 




M b 


ref c 






[M ] 


[log(M Q yr- 1 )] 




f^V TV. 

Y lau 


1V11 


U.OU 


-o. id 


1 

1 


DE Tau 


A/TO 

IVlz 


U.OO 


- / .y / 


1 

1 


lau 


1V1U 


U. /U 


- / . lo 


1 


vjl lau 


IV u 


U.oO 


-o.Uo 


1 

1 


r^A/f Aiiv 
^ivi Am 


TY"7 


n 7Q 


Q 11 

-o. 11 


1 

1 


1"P Tan 

ijr lau 


A/rn 

1V1U 


n 7n 
U. i U 


-y.uo 


1 

1 


uu lau A 


A/TQ 

IVlo 


n Q7 
U.o / 


-y. i 


n 

z 


DF Tau 


M0.5 


0.68 


-8.0 


2 


FO Tau A 


M2 


0.60 


-7.9 


2 


T?V Tern A 

r v lau I\ 


IV 


i in 
1 . lz 


-(.( 


o 

A 


C*C* Ton A q 

Li v_j iau /\a 


l\ I 


n 7K. 

U. i 


7 Q 

- / .y 


n 
A 


— Ab 


M0.5 


0.68 


-8.2 


2 


— Ba 


M6.0 


0.12 


-10.5 


2 


GH Tau A 


M1.5 


0.63 


-8.9 


2 


- B 


M2 


0.56 


-9.4 


2 


IS Tau A 


K7 


0.79 


-8.1 


2 


RW Aur A 


Kl 


1.34 


-7.5 


2 


- B 


K5 


0.93 


-8.8 


2 


T Tau A 


KO 


2.11 


-7.5 


2 


UZ Tau A 


Ml 


0.65 


-8.7 


2 


— Ba 


M2 


0.56 


-8.8 


2 


V807 Tau A 


K7 


0.76 


-8.4 


2 


V955 Tau A 


K5 


0.87 


-8.5 


2 


XZ Tau A 


M3 


0.44 


-8.9 


2 


CIDA-14 


M5 


0.15 


-10.3 


3 


MHO-6 


M4.75 


0.18 


-10.3 


3 


V410 Anon 13 


M5.75 


0.11 


-11.3 


3 



a Masses from Lyon tracks. 

b Accretion rates from veiling, UV-excess and Ha profile mod- 
cling, in references [1], [2] and [3] respectively (see note [c] below). 

c References for the spectral types and M . 1: Gullbring et al. 
1998; 2: WG01; 3: MHCBH03. The Gullbring ct al. M , calcu- 
lated from t/-band veiling using DM tracks, have been modified 
here to Lyon masses and radii; the WG01 M , from UV-excess 
but using Lyon tracks, have been adopted unchanged. The 
MHCBH03 M from Ho modeling are also adopted unchanged, 
but the masses given in in MHCBH03, from DM tracks, have 
been changed to Lyon ones. See §5.3. 



